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ABSTRACT 


The  application  of  LORAN-C  in  the  hyperbolic  mode  as  a 
positioning  systes  for  hydrographic  surveys  was  investi¬ 
gated.  Observed  LORAN-C  tine  differences  froa  a  field  test 
conducted  in  Monterey  Bay,  California  were  compared  to 
calculated  tiae  differences  determined  from  geographic  posi¬ 
tions  based  on  a  microwave  positioning  system.  Four  methods 
were  used  to  determine  the  calculated  tiae  differences.  The 
first  three  methods  were  (1)  applying  only  the  seawater 
Secondary  Factor,  (2)  computing  the  time  difference  based  on 
a  Semi-Empirical  TD  Grid,  and  (3)  applying  ASF  Correctors 
from  the  DMAHTC  LORAN-C  Correction  Table.  The  final  method 
applied  multiple  observed  ASF  Correctors  at  five  minute 
latitude  and  longitude  intervals.  By  applying  multiple 
observed  ASF  Correctors,  which  was  the  most  accurate  method, 
a  38.3  merer  1  drms  with  a  lane  offset  of  3  to  12  meters 
using  the  9940  X-T  LORAN-C  combination  was  obtained.  Based 
upon  the  results  presented,  it  may  be  possible  to  use 
LORAN-C  for  hydrographic  surveys  at  scales  of  1:  80,000. 


TABLE  0?  CONTENTS 


a 

I.  I  NT  HOD  OCT  ION . 9 

A.  OSE  OP  LOBAN-C . 9 

.  B.  APPLICATION  OP  ASP  CORRECTORS  POR  NAVIGATION  .  11 

|  '  C.  APPLICATION  OF  ASF  TO  HYDROGRAPHIC  POSITIONING  13 

D.  OBJECTIVES . 14 

II.  NATURE  OP  THE  PROBLEM . 16 

A.  THE  PRINCIPLES  OP  LORAN-C . 16 

|  B.  PHASE  LAG . 19 

C.  TD  MODEL . 22 

D.  SEMI -EMPIRICAL  TD  GRID  CALIBRATION  MODEL  ...  24 

1.  Technical  Approach . 25 

i  2.  Generalized  Range/Bearing  Model  .  27 

E.  DMAHTC  MODEL . 31 

1.  Sea  SF  Modal . 31 

2.  Land  SP  Model . 31 

|  3.  Table  Description . 34 

F.  ATTEMPTED  DETERMINATION  OP  ASF  CORRECTORS  BY 

HYDRO  FIELD  PARTIES . 38 

III.  EXPERIMENTAL  PROCEDORE  .  42 

)  A.  ..ELD  PROCEDURES . 42 

B.  MICROWAVE  SYSTEM  POSITIONING  .  45 

C.  LORAN-C  COMPUTATIONS . 46 

1.  Seawater  Secondary  Factors  (SF) . 47 

?  2.  Semi-Empirical  TD  Grid . 49 

3.  Calculated  Table  ASF  Correctors  .  51 

4.  Observed  ASF  Correctors . 52 

IV.  RESULTS . 54 

1  A.  SEAWATER  SECONDARY  FACTORS  .  55 

B.  SEMI -EMPIRICAL  T  D  GR  ID  . . 57 

C.  TABLE  ASF  CORRECTORS . 59 

D.  MULTIPLE  OBSERVED  CORRECTORS . 60 

i 


5 


V 


CONCLUSIONS 


APPENDIX  A:  PROGRAM  GPBYLQ 

APPENDIX  B:  PROGRAM  LOR  AN 

APPENDIX  C:  PROGRAM  LOPLC 

APPENDIX  D:  PROGRAM  LORTAB 

APPENDIX  E:  PROGRAM  ASFSEL 

BIBLIOGRAPHY  . 

INITIAL  DISTRIBUTION  LIST  . 


LIST  OF  TABLES 


Phase  Retardation  or  Lag  of  Radio  Raves 
Coding  Delay  9940  LORAN-C  Chain  .  .  .  . 

TD  Bias  (b)  -  psec  . 

Seawater  Coefficients . 

Land  Coefficients . 

Sea  SF  Model  Coefficients . 

Geographic  Names  and  Positions . 

Microwave  Positioning  Equipment  .  .  .  . 

LORAN-C  9940  Chain  Data . 

Data  Set  Parameter  . 

Seawater  Secondary  Factor  Errors  .  .  .  . 
Semi-Empirical  TD  Grid  Correction  Errors 
Multiple  LORAN-C  Correction  Table  Errors 
Multiple  Observed  Correction  Errors  .  . 


LIST  OF  FIGURES 


Hyperbolic  Fix  (From  Coast  Guard  LORAH-C  User 

Handbook,  1974)  . . 

Transmitted  Radiation . 

Phase  Lag . . 

Location  of  West  Coast  LORAN-C  Stations  .  .  .  . 

Hixed  Path  TD  Geometry . 

Composite  Land-Sea  Path  . 

99U0-W  ASF  Correction  Table . 

9940-X  ASF  Correction  Table  . 

9940-Y  ASF  Correction  Table  . 

Location  of  Test  Areas  in  Monterey  Bay  . 


I.  INTRODUCTION 


A.  USE  OF  LORAN-C 

In  recent  years  there  has  been  an  increasing  usage  of  a 
LORAN-C  receiver  and  the  LORAN-C  network  as  the  primary 
horizontal  control  for  such  scientific  studies  as  deep  ocean 
dumpsites,  marine  fisheries  studies  [Rulon,  1979],  bathyme¬ 
tric  surveys,  and  recently,  a  reconnaissance  hydrographic 
survey.  Examples  of  bathymetric  surveys  conducted  by  the 
National  Oceanic  and  Atmospheric  Administration  (NOAA) , 
which  have  used  Loran-C  for  positioning  are: 

1)  Su- 100-1-79  Gulf  of  Alaska  [NOAA  H-9822,  1979], 

2)  SU- 100-2-79  Gulf  of  Alaska  [NOAA  H-9823,  1979], 

3)  S-D902-WH-82  U.S.  West  Coast  [NOAA  Ship  Surveyor, 
1982]. 

Recently  an  attempt  was  made  to  use  LORAN-C  as  the 
sounding  position  control  for  a  reconnaissance  hydrographic 
survey  S-K902-WH-82  [NOAA,  1982  ].  This  was  a  special  survey 
conducted  by  the  NOAA  Ship  Whiting  in  May  1982.  Special 
surveys  are  field  examinations  of  very  limited  extent  or 
scope  and  frequently  require  unique  survey  or  data  collec¬ 
tion  procedures  [Umbach,  1976  ].  The  purpose  of  this  project 
was  to  verify  the  existence  and  extent  of  reported  shoaling 
in  three  safety  fairways  in  the  Gulf  of  Mexico  [NOAA,  1982], 

The  use  of  Loran-C  as  a  positioning  system  for  basic 
hydrographic  surveys  has  been  very  limited  due  to  the  abso¬ 
lute  accuracy  of  the  long  range  system.  A  basic  hydro- 
graphic  survey  is  defined  as  a  survey  which  is  so  complete 
that  it  need  not  be  supplemented  by  other  surveys.  "It  must 
be  adequate  to  supersede  for  charting  purposes  all  prior 


surveys"  [Umbach,  1976].  Variables  which  affect  the  accu¬ 
racy  of  LORAN-C  are  signal  propagation  variations,  weather, 
and  sky  waves.  The  effects  of  weather  and  sky  waves  on 
LORAN-C  propagation  are  best  described  by  Samaddar  [1980] 
and  the  American  Practical  Navigator,  [DMA,  1977]  respec¬ 
tively. 

Signal  propagation  variations  are  due  to  the  phase 
retardation  of  the  signal  as  it  passes  over  an  all  se?  ater 

path,  over  land  paths,  or  partial  seawater-land  pa  s  as 

compared  to  free  space.  Table  I  summarizes  phase  re  _ la- 

tion  changes  [Mortimer,  1  978],  Errors  due  to  an  a  sea 


TABLE  I 

Phase  Retardation  or  Lag  of  Radio  Waves 


Propagation  Path 

Sep  res  en  tat  ion 
Propagation 
Velocity 
(km /sec) 

Difference  in  Phase 
Lag  at  500  km 
Compared  with 

Wave  in  Line 

Above  (m) 

Vacuum 

299792.5 

- 

Direct  wave  through 
earth's  atmosphere 

299691 

170 

Ground  wave  over 
sea  water 

299560 

220 

Ground  wave  over 
rugged  mountains 

298  899 

1,300 

water  path  are  known  as  the  Secondary  Factor  (SF)  and  errors 
due  to  a  land  path  or  aired  path  are  known  as  the  Additional 
Secondary  Factor  (ASF)  [Speight,  1982  ]. 


ASF  Corrections  in  tha  LORAN-C  system  can  be  as  large  as 
pins  or  minus  four  microseconds,  which  is  600  meters  on  the 
baseline.  In  other  areas  with  the  same  LORAN-C  coverage, 
these  corrections  aay  be  much  larger  due  to  the  expansion  of 
the  distance  between  adjacent  hyperbolic  lines  of  position. 
For  example,  at  32®  N  and  80®  W ,  using  lattice  pair  9960-X, 
a  four  microsecond  ^isec)  error  will  offset  the  9960-X  line 
of  position  approximately  2438  meters  [Speight,  1982  ]. 

B.  APPLICATION  OF  ASF  CORRECTORS  FOR  NAVIGATION 

To  compensate  for  the  Loran-C  positional  errors  caused 
by  the  ASF  Correctors,  the  Secretary  of  Transportation 
tasked  the  Defense  Mapping  Agency  Hydrographic/Topographic 
Center  (DMAHTC)  ,  the  National  Oceanic  and  Atmospheric 
Administration  (NOAA)  ,  and  the  United  States  Coast  Guard 

(USCG)  with  the  job  of  determining  and  applying  the 

Additional  Secondary  Correctors  for  each  Loran-C  chain. 
This  task  was  published  in  the  Department  of  Transportation 
(DOT)  National  Plan  for  Navigation  in  the  July  19,  1974 

Federal  Register.  These  corrections  should  provide  95% 

assurance  that  a  vessel  could  fix  its  position  to  a  pred¬ 
icted  accuracy  of  1/4  nautical  mile  (NM)  within  the  U.  S. 
Coastal  Confluence  Zone  (CCZ)  and  the  Great  Lakes.  The  CCZ 
is  defined  as: 

"the  seaward  approaches  to  land,  tve  inner  boundary  of 
which  is  the  harbor  entrance  and  the  outer  boundary  of 
which  is  50  nautical  miles  offshore  or  the  edge  of  the 
Continental  Shelf  (100  fathom  contour)  whichever  is 
greater.  " 

The  1/4-NM  accuracy  requirement  also  affects  the 

nautical  chart.  The  National  Ocean  Survey  (NOS)  ,  which 
publishes  charts  for  the  CCZ,  engaged  in  a  program  with  the 
USCG  and  DMAHTC  to  provide  the  coastal  navigator  with  charts 
overprinted  with  lattices  which  meet  1/4-NM  accuracy.  The 
USCG,  as  operator  cf  the  LORAN-C  radior.avigation  system. 


conducts  surveys  to  ensure  that  LORAN-C  coverage  exists 
within  the  CCZ  and  will  be  reponsible  for  the  verification 
of  1/4-NM  accuracy  for  all  coastal  LORAN-C  service.  In 
conjunction  with  NOS,  it  assists  in  surveys  of  coastal 
waters  of  the  United  States  to  allow  production  of  LORAN-C 
charts  based  on  observed  field  data  to  aeet  the  standards 
set  forth  above  (Speight,  1  982]. 

DHAHTC ,  for  LORAN-C  civil  need,  prepares  grid  predic¬ 
tions  froa  its  data  base.  Based  oa  analysis  and  verifica¬ 
tion  of  the  predicted  grid  from  a  USCG  and/or  NOS  survey,  it 
produces  revisions  to  the  initial  grid  predictions  [Speight, 
198  2  ].  At  present,  DHAHTC  has  provided  NOS  with  ASF 
Corrected  LORAN-C  Lattices  which  are  overprinted  on  the  NOS 
Charts.  Each  chart  with  ASF  Correctors  applied  contains  one 
of  the  following  notes: 

"The  LORAN-C  lines  of  Dosition  overprinted  on  this  chart 
have  been  prepared  for  use  with  groudwave  signals  and 
are  presently  compensated  only  for  theoretical  propaga¬ 
tion  delays,  which  have  not  yet  been  verified  by 
observed  data.  Hariners  are  cautioned  not  to  rely 
entirely  on  the  lattices  in  inshore  waters.  Skywave 
corrections  are  not  provided". 


or 


"The  LORAN-C  lines  of  position  overprinted  on  this  chart 
have  been  prepared  for  use  with  groundwave  signals  and 
are  compensated  with  propagation  delays  computed  from 
observed  data.  Mariners  are  cautioned  not  to  rely 
entirely  on  the  lattices  in  inshore  waters.  Skywave 
corrections  are  not  provided"  [Speight,  1982]. 

Presently,  all  of  the  NOS  Charts  of  1:80,000  to 
1:120,000  scale  covering  the  east  coast.  Gulf  coast,  and 
Great  Lakes  show  LORAN-C  lattices  that  have  been  compensated 
for  Additional  Secondary  Factors.  Most  of  the  lattices  on 
these  charts  have  been  constructed  from  DMAHTC  data  tapes 
that  provide  adjusted  LORAN-C  readings  for  each  rate  at 
every  five  minutes  of  latitude  and  longitude.  A  few 
lattices  were  constructed  using  a  single  ASF  Correction  for 
the  entire  chart  area.  Five  minute  data  tapes  were  not 
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furnished  by  DMAHTC  for  constructing  lines  of  position  for 
LORAN-C  rates  on  the  Best  C  xst  Charts.  On  these  charts  a 
single  average  ASF  Correction  was  used  to  adjust  each 
lattice  [NOAA,  Marine  Chart  Division,  1982]. 

In  addition  to  supplying  corrected  LORAN-C  lattices  for 
nautical  charts,  DMAHTC  prepares,  distributes,  and  periodi¬ 
cally  updates  unclassified  ASF  LORAN-C  Correction  Tables 
[Speight,  1982].  The  ASF  Correction  Tables  are  for  preci¬ 
sion  navigation,  utilizing  digital  computers  to  convert 
LORAN-C  time  differences  to  geographic  coordinates  [Speight, 
1982].  Presently,  the  ASF  correctors  found  in  the  LORAN-C 
Correction  Tables  were  determined  using  theoretical  propaga¬ 
tion  delays.  ASF  correctors  listed  in  the  tables  are  going 
to  be  updated  with  observed  data  and  reprinted  the  first 
quarter  of  1983  [Rallace,  1  982]. 

C.  APPLICATION  OF  ASF  TO  HYDROGRAPHIC  POSITIONING 

Schr.ebele  [  1979]  investigated  the  possibility  of  using 
Loran-C  as  an  electronic  positioning  system  for  hydrographic 
surveying.  He  concluded  that  in  Monterey  Bay,  California  a 
single  Additional  Secondary  Factor  (ASF)  applied  to  offshore 
lines  of  position  gave  a  root  mean  square  error  (drms)  of  66 
meters  for  the  Best  Coast  9940  Y-B  pair  and  a  predicted  42 
meter  drms  error  for  9940  X-Y  rates. 

The  42  meter  predicted  drms  is  larger  than  Nelson's 
[General  Electric  Co.,  197  9]  findings  in  San  Francisco  Bay. 
He  demonstrated,  in  a  dynamic  mode,  that  the  precision  of 
LORAN-C  was  60.8  meters  2  drms  (30.4  meter  1  drms)  with  a 
worst  case  of  71.2  meters  2  drms  (35.6  meter  1  drms).  A 
mean  difference  or  offset  between  the  measured  time  differ¬ 
ence  and  the  calculated  time  difference  was  34  nanoseconds 
for  the  9940-X  rate  and  one  nanosecond  for  the  9940-Y  rate. 
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He  also  obtained  a  precision  of  38.0  aeters  2  drms  (19.0 
Deters  1  drms)  in  the  static  mode.  Nelson  also  states,  that 
the  above  precision  is  only  achievable  if  the  user  has  a 
LORAN-C  receiver  which  has  the  performance  capabilities  of 
those  used  in  the  experiment.  The  LORAN-C  receiver  must 
have  "comparable  signal  averaging  rime,  extra  notch  filters, 
and  attenuation  of  the  signal"  [General  Electric  Co.,  1979]. 

D.  OBJECTIVES 

The  National  Ocean  Survey  requires  that  hyperbolic 
control  systems  used  for  hydrographic  surveying  exhibit  a  1 
drms  of  less  than  0.5  millimeter  at  scale  of  the  survey 
[Ombach,  1976  ].  Although  this  requirement  is  generally  for 
2  mHz  phase  comparision  systems,  it  can  be  inferred  that  if 
also  applies  to  other  hyperblic  systems  such  as  LORAN-C. 
The  scale  routinely  used  for  coastal  surveys  is  between 
1:40,000  and  1:80,000  [Ombach,  1976  ]  yielding  an  allowable 
error  of  20  to  40  meters  not  including  systematic  errors. 
Schnebele  [1979]  concluded  that  hyperbolic  LORAN-c,  after 
applying  a  single  ASF  Corrector,  is  unsuitable  for  basic 
hydrographic  surveying. 

Whether  or  not  applying  multiple  Additional  Secondary 
Factors  (ASF)  to  LORAN-C  lines  of  position  will  reduce  the 
drms  sufficiently  to  meet  the  accuracy  standards  set  by  the 
National  Ocean  Survey  Hydrographic  Manual  will  be  ascer¬ 
tained  in  this  study.  The  term  multiple  ASF  Correctors 
refers  to  the  application  of  more  than  one  corrector  to 
LORAN-C  lines  of  position  over  a  given  area.  The  variable 
ASF  Correctors  result  from  varying  delays  of  the  electromag¬ 
netic  wave  as  it  propagates  over  different  land  segments. 

Three  methods  cf  applying  multiple  ASF  Correctors  were 
tested.  The  first  method  was  the  application  of  a 


Sea i-Eapir ical  Tiae  Difference  Grid  Calibration  Model  devel¬ 
oped  by  The  Analytic  Science  Corporation  [  1979].  The  spon¬ 
soring  agency  vas  the  United  States  Coast  Guard.  The 
Sea i-Eapir ical  Model  applies  Secondary  Factors  and 
Additional  Secondary  Factors  for  each  geodetic  position 
based  on  the  distance  over  land,  the  distance  over  water, 
and  the  total  distance  using  aean  sea  water  and  land  conduc¬ 
tivities. 

The  second  aethod  which  was  investigated  applies  ASF 
Correctors  found  in  the  DHAHTC  LORAN-C  Correction  Tables 
[DMAHTC,  1981]  to  LORAN-C  lines  of  positions.  These  correc¬ 
tors  were  derived  froa  the  ground  conductivities  which  have 
been  deterained  in  the  field  by  a  Coast  Guard  calibration 
team  [tJ.S.  Naval  Oceanographic  Office,  1982]. 

Finally,  a  third  aethod  was  pursued.  ASF  Correctors, 
which  were  determined  by  field  observations,  were  applied  to 
the  LORAN-C  lines  of  position.  These  were  determined  by 
computing  the  difference  between  the  observed  LORAN-C  rates 
and  the  expected  time  difference  which  was  calculated  using 
four  lines  of  position  from  a  very  accurate  microwave  posi¬ 
tioning  system.  These  ASF  Correctors  were  determined  at 
five  minute  latitude  and  longitude  intervals. 


zi.  s i  xu  2i asi5fi 

A.  THE  PRINCIPLES  OP  LO R1H-C 

To  understand  the  problems  associated  with  LORAN-C  when 
used  during  hydrographic  operations,  one  oust  first  under* 
stand  its  principles  of  operation.  LORAN-C  is  a  low 
frequency,  pulsed  signal,  hyperbolic,  radio  navigation 
system,  eaploying  tiae  difference  aeasureaents  of  signals 
received  by  the  navigator  from  at  least  three  ground  tran¬ 
smitting  stations  [Speight,  1982].  The  stations  are 
comprised  of  a  master  transmitting  station,  two  or  more 
secondary  transmitting  stations  which  are  strategically 
spaced  several  hundred  miles  apart  and,  if  necessary,  a 
System  Area  Noniter  (SAN)  Station  [O.S.  Coast  Guard,  1974], 

System  Area  Honitor  (SAH)  stations  associated  with  each 
LORAH-C  chain  apply  differential-type  corrections  to  the 
rates  in  real-time.  SAN  stations  continuously  monitor  the 
signals  from  all  transmitters  in  the  chain.  If  the  oberved 
tiae  difference  deviates  by  -  more  than  0.05  jusec  from  the 
expected  value,  then  the  appropriate  secondary  adjusts  its 
emission  delay  time  in  order  to  remove  the  error  [Schnebele, 
1979]. 

The  master  and  at  least  two  secondary  stations  are 
located  such  that  the  signals  from  the  transmitting  stations 
can  be  received  throughout  the  desired  coverage  area.  The 
master  station  is  designated  by  the  letter  "H"  and  the 
secondary  stations  or  slave  stations  are  desigr.atied  W,  X, 
I,  or  Z  [.O.S.  Coast  Guard,  1974  ]. 

All  stations  transmit  on  the  common  frequency  of  100 
kHz.  Interference  between  transmitters  is  avoided  through 
the  use  of  time  separation  [Poppe,  1982  ].  After  the  master 
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station  transmits  a  pulse,  each  secondary  station  delays  its 
own  transaission  for  a  fixed  tine,  called  the  secondary 
coding  delay.  This  coding  delay  is  synchronized  through  the 
use  of  cesiua  frequency  standards  at  each  station.  The  high 
stability  and  accuracy  of  these  standards  perait  each 
station  to  derive  its  ovn  tiae  of  transaission  without 
reference  to  another  station  [ DMA,  1977  ],  Secondary  coding 
delays  are  predeterained  by  system  propagation  times  and 
equipment  characteristics  [  Laurila,  1976  ]. 

The  pulse  from  the  master  transmitter  is  distinguished 
from  those  of  the  secondaries  through  phase  coding  of  the 
pulses.  Phase  coding  refers  to  the  inversion  of  the  nega¬ 
tive  and  positive  peaks  of  the  sine  wave  comprising  the  100 
kHz  carrier  portion  of  the  pulse.  The  purpose  of  the  phase 
coding  is  twofold: 

"First,  it  permits  automatic  discrimination  between  the 
master  and  the  various  secondary  stations,  thereby 

fermittmg  all  stations  to  be  identified  by  their  rela- 
ive  timing  with  respect  to  the  master" 

"Second,  the  phase  coding  provides  protection  against 
excessively  long  skywave  delays  which  would  cause  the 
late  arrival  of  the  preceeding  pulse  to  coincide  with 
the  leading  edge  or  groundwave  portion  of  a  pulse  being 
tracked"  [Poppe,  1982]. 


The  signals  are  received  by  a  mobile  receiver  where  the 
differences  in  tiae  of  arrival  of  the  master  signal  and 
various  secondary  signals  are  measured  and  displayed  on  the 
indicator  portion  of  the  LORAN-C  set.  The  accuracy  of  this 
tiae  difference  is  increased  by  phase  comparision  "of  the 
synchronized  100  kHz  carrier  within  the  master  and  secondary 
pulses"  [Laurila,  1976].  This  measured  tiae  difference  (TD 
-  in  microseconds)  represents  a  hyperbolic  Line  of  Position 
(LOP)  [U.S.  Coast  Suard,  1974],  The  intersection  of  two  or 
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Hyperbolic  Pix  (From 
Handbook,  1974) 


Coast  Guard  LORAN-C  User 


more  LORAN-C  LOP’s  defines  the  position  of  the  observer 
(Figure  2.1).  When  plotted  on  a  chart,  the  intersection  of 
the  resultant  hyperbolic  lines  defines  a  geographical 
position  f  Speight,  1982]. 


B.  PHASE  LAG 


In  a  vacuum,  the  velocity  of  radiated  energy  from  an  antenna 
for  LORAN-C  is  299792.458  ka/sec.  Since  radiated  energy 
cannot  be  shown  pictorially,  the  phase  of  the  transmitted 
radiations  is  used.  The  lines  of  constant  phase  of  the 
transmitted  radiation  are  shown  in  Pigure  2.2  by  the  curved 
lines  labeled  aa',  bb' ,  and  other  similar  designations. 
They  define  the  wave  front  as  it  proceeds  outward  froa  the 
antenna  in  all  directions.  The  distance  between  each  line 
of  constant  phase  is  one  wavelength  (x)  [Admiralty  Manual, 
196  5]. 

X  (meters)  *  29  979  2.4  58  km/sec  *  frequency  in  kHz 


Antenna 


Figure  2.2  Transmitted  Radiation 


The  velocity  of  the  radiated  energy  in  air  depends  on 
eaperature,  pressure,  humidity,  and  the  nature  of  the 
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surface  over  which  the  transmissions  pass.  The  retardation 
of  a  transmitted  wave  is  known  as  phase  lag.  When  low 
frequencies  are  employed,  such  as  L03AN-C  at  100  kHz,  the 
effects  of  change  of  temperature,  pressure,  and  huaidity  are 
swaaped  by  the  effects  caused  by  changes  in  the  nature  of 
the  surface  over  which  the  transaissions  are  traveling 
[Admiralty  Manual,  1965].  The  General  Electric  Company, 
TEMPO  division,  conducted  a  LORAN-C  Signal  Analysis 
Experiment  under  the  direction  of  the  O.S.  Coast  Guard. 
This  experiment  was  conducted  along  tht*  O.S.  West  Coast. 
The  General  Electric  Company  recorded  an  overall  change  of 
108  nsec  and  116  nsec  time  of  arrival  from  the  master  and 
X-secondary  stations  respectively  after  a  storm  [Samaddar, 
1980].  If  ASF  Correctors  are  as  large  as  two  microseconds 
on  the  West  Coast  [DMAHTC,  1981]  then  the  ASP  Corrections 
are  20  times  larger  in  magnitude  than  weather  effects  for 
the  9940  West  Coast  LORAN-C  chain. 

The  change  in  transmission  rates  or  phase  lag  are  a 
result  of  the  amount  of  energy  transferred  from  •‘•he  tran¬ 
smitted  radiation.  This  energy  transfer  depends  on  the 
absorption  qualities  (inversely  related  to  conductivity)  of 
the  surface  over  which  they  are  passing  and  their  wavelength 
(or  frequency).  The  lower  the  conductivity  and  the  longer 
the  wavelength  (or  lower  the  frequency)  the  greater  the 
transfer  of  energy,  and  vice  versa.  Seawater  has  a  rela¬ 
tively  high  conductivity  (5.0  mhos/meter) .  Land  has  a  much 
lower  conductivity,  which  varies  from  marsh  (fairly  high)  to 
dry  sand  and  rock  (very  low)  [Admiralty  Manual,  1965  ].  Two 
excellant  papers  that  discuss  the  electrical  properties  of 
soil  are  those  of  Smith-Rose  [  1934  ]  and  Pressey,  Ashwell, 
and  Fowler  [  1956].  Smith-Rose  [  1934  ]  found  that  the  conduc¬ 
tivities  for  soil  ranged  from  0.18  mhos/meter  for  a  grey 
clay  with  salt  to  .00001  mhos/meter  for  granite. 


Sea  Level 


Distance  from  Antenna  (k as) 


Figure  2.3  Phase  Lag 


Phase  lag  is  illustrated  in  Figure  2.3.  The  lines  of 
constant  phase  ,  aa',  fab',  and  cc*  become  distorted  as  they 
progress  along  the  sea  surface.  The  dotted  lines  represent 
the  position  of  the  lines  of  constant  phase  in  the  absence 
of  the  sea  surface.  The  wavelength  (X)  in  meters,  measured 
at  heights  of  several  wavelengths  above  the  sea,  remains 
about  the  same  as  the  direct  wave  through  the  earth's  atmo¬ 
sphere  at  299691  km/sec  *  frequency  in  kHz.  Near  sea  level 
the  absorption  of  energy  retards  the  progress  of  the  wave- 
front,  and  makes  the  sea  level  wavelength  (x‘)  less  than  A. 
As  the  lines  of  constant  phase  progress  away  from,  the 
antenna  the  phase  lag  increases  with  distance.  This  is 
known  as  the  Secondary  Factor  (SF)  [Admiralty  Manual,  1965]. 

The  most  intriguing  feature  of  phase  lag  occurs  at  the 
coastal  interface  where  there  is  an  extreme  change  in  the 
conduct ivit es  between  land  and  sea.  Visualizing  the  wave 
front  in  three  dimensions,  the  lower  part  of  the  wave, 
slowed  by  the  drag  of  the  ground,  lags  further  and  further 
behind  the  upper  part  as  the  wave  crosses  the  land.  At  the 


coastline  it  suddenly  encounters  the  much  lower  impedance  of 
the  sea,  and  in  a  very  short  distance  the  bottom  of  the  wave 
tries  to  catch  up  with  the  top,  as  though  the  whole  wave 
front  were  an  elastic  balloon.  This  is  known  as  "phase 
recovery"  [Eaton,  1979].  Phase  recovery  was  verified  during 
tests  on  Decca  transmissions  across  the  south  coast  of 
England  by  Pressey,  Ashwell,  and  Fowler  [  1956]. 

The  determination  of  the  Secondary  Factor  for  seawater 
is  fairly  direct  since  the  conductivity  of  seawater  (5.0 
mhos/meter)  is  fairly  constant.  But  for  land  the  conduc¬ 
tivity  can  vary  depending  on  the  type  of  soil  and  its  water 
content  [Smith-Rose,  1934J. 

Phase  lag  for  radiated  energy  over  land  can  be  deter¬ 
mined  two  ways: 

1)  Assign  an  average  land  conductivity  to  the  ASF 
Model.  For  example,  the  average  conductivity  for 
the  soil  on  the  west  coast  is  0.003  mhos/meter. 

The  average  land  conductivity  will  determine  the 
average  phase  retardation  of  the  path  [The 
Analytic  Science  Corporation  ,  1979  ]. 

2)  Determine  every  conductivity  for  each  portion 
of  a  line  segment  from  the  transmitter  to  the 
receiver.  The  total  of  these  conductivity 
segments  constituting  a  land-water  profile  will 
determine  the  total  phase  retardation  cf  the  parh 
[Speight ,  1982  ]. 

C.  TC  MODEL 

Positional  fix  accuracy  using  LOHAN-C  is  primarily 
dependent  on  a  chart  makers  ability  to  accurately  compute 
the  expected  difference  in  time-of-arrival  (TOA)  of  received 
groundwave  signals  from  the  transmitting  s+ations.  Time 


differences  (TD) ,  are  the  differences  between  the  TOAs  of 

the  secondary  and  saster  transmitters. 

TD.  »  TOA  -  TOA  (2.1) 

lie 

i  *  secondary  Station 
a  =  Has  ter  Station 


TOA  coapatations  are  dependent  upon  an  accurate  knowledge  of 
the  signal  phase  delay. 


The  phase  delay 

of  a  groundwave 

signal  is  generally 

expressed  as: 

♦  a 

T  ♦  SF 

(2.2) 

nR 

3 

♦  SF 

C 

where  n  is  the  surface  refractive  index,  C  is  the  speed  of 
light  in  a  vacuum,  H  is  the  range  between  the  receiver  posi¬ 
tion  and  the  transmitting  station.  The  primary  phase  delay, 
T,  is  the  computed  travel  time  of  the  LORAH-C  pulse  over  a 
distance  equal  to  the  tran saitter-to-receiver  great  circle 
path  length,  taking  into  account  the  velocity  of  electromag¬ 
netic  waves  and  the  index  of  refraction  of  the  atmosphere. 
The  secondary  factor  (SF)  is  a  correction  to  the  primary 
phase  delay  and  accounts  for  the  phase  lag.  The  dominant 
term  in  (2.2)  is  the  primary  phase  delay  (T) .  The  SF  is 
usually  an  order  of  magnitude  smaller  [The  Analytic  Science 
Corporation,  1979]. 

Thus,  time-of -arrivals  can  be  expressed  as: 
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where  CD  is  the  true  emission  delay  or  coding  delay  for  the 
LORAN-C  chain  [The  Analytic  Science  Corporation,  1979].  The 
coding  delay  is  equal  to  a  time  delay  plus  a  computed  one 


way  baseline  tine  (Be)  which  includes  the  secondary  phase 
correction  for  an  all  seawater  path.  The  oneway  baseline 
time  (Be)  is  equal  to  the  distance  between  the  master  and 
secondary  transmitters  in  neters  divided  by  the  propagation 
velocity  of  LORAN-C  through  the  earth* s  atmosphere  (299.691 
meters  per  microsecond  [  Navigation  Department  DMA,  1982]). 
See  Table  II  for  Coding  Delay  values  for  the  9940  chain 

TABLE  II 

Coding  Delay  9940  LORAN-C  Chain 


I 


[Riordan,  1979].  Combining  equations  2.1,  2.3,  and  2.4,  the 
true  TD  is  given  by  equation  2.5  [The  Analytic  Science 
Corporation,  1979]. 

TD  =  (T  -  T  )  ♦  (SF  -  SF  )  +  CD  (2.5) 

i  i  m  i  m  i 

D.  SEMI -EMPIRICAL  TD  GRID  CALIBRATION  MODEL 

The  Semi-Empirical  TD  Model  was  developed  by  The 
Analytic  Science  Corporation  [  1979]  in  Reading, 
Massachusetts  for  the  West  Coast  9940  LORAN-C  chain. 
Similar  "time  difference  (TD)  grid  calibration  techniques 
have  been  successfully  employed  to  develop  an  accurate 
(approximately  100  nsec  drms)  calibrated  grid  for  St.  Marys 
River  LORAN-C  chain",  [The  Analytic  Science  Corporation, 


CD  ♦  Be  =  1  1000  ♦  2796.90  =  13796.90  jisec 
CD  ♦  Be  =  2  7000  +  1094.49  =  28094.49  jisec 
CD  +  Be  =  4  0000  ♦  1967.27  =  41967,27  /isec 


Pair  9940-W 
994 0-X 
9940-Y 


1 •  Technical  App roach 


The  Semi-Empirical  Model  is  based  cn  Millington's 
empirical  approach  for  computing  the  secondary  factor  over  a 
mixed  (multiple-homogeneous  segment)  path  which  combines 
land  and  sea  phase  delays.  The  generalized  semi-empirical 
polynomial  functional  form  for  the  SP  of  the  LOSAN-C  station 
is  given  by: 


SF 


j 


where 


K 

2  L 

SF  (T  ,  6  )  *  Z  AT  +  S(C  Sial  S  +  D  cosl  0  )  (2.6) 

j  j  k=-Ki  k  j  1=1  jl  j  jl  j 


j  =  secondary  (W,X,  or  Y)  or  master  (M)  station, 
nB 


T  = 

j 

R  = 

j 

8  = 
j 


C 

jl 


-  =  jth  station-to  user  primary  phase  delay 

c 

jth  station -to -user  great-circle  path  length, 
user  path  bearing  angle  at  the  jth  station, 
and  L  are  positive  integers, 
and  D  are  the  station-dependent  coefficients 

jl 

of  harmonic  terms  in  the  model, 
is  the  range-dependent  coefficient  of  the  model 

which  may  in  general  be  station-dependent. 
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Data  from  27  coastal  sites  distributed  along  the 

i 

West  Coast  and  122  land-sea  sites  distributed  in  the 
Southern  California  CCZ  (between  Point  Arguallo  and  San 
Diego  -  see  Figure  2.4)  were  used  in  a  Kalman  estimation 
algorithm  to  compute  the  uncertain  coefficients  of  the  land 
and  sea  models  of  the  ID  grid  calibration  algorithm.  (An 
explanation  of  Kalman  filtering  for  the  layman  is  presented 
by  Roger  M.  du  Plessis  [  1  967].)  The  calibrated  algorithm 


vas  used  to  compute  TDs  at  each  data  site  and  the  TD  resi¬ 
duals  (difference  between  measured  and  calibrated  TDs)  were 
examined.  Adjustments  were  then  made  to  the  TD  model  struc¬ 
ture  in  an  attempt  to  further  reduce  the  residuals.  This 
process  of  adjusting  the  model  structure  is  repeated  until 
the  residuals  agree  with  the  expected  theoretical  covariance 
associated  with  the  TD  model.  The  model  which  exhibited  the 
'•best"  performance  was  selected  as  the  West  Coast  TD  grid 
model. 

2.  Generalized  R ana e/B earing  Model 

The  Generalized  Range/Bearing  (GRB)  Model  was 
selected  as  the  "best11  semi-empirical  calibration  model  for 
the  West  Coast  chain.  The  semi-empirical  function  is: 

L 

SF  =  A 0  *  AIT  +  e  (C  sinlB  ♦  D  cosl  3  )  nsec  (2.7) 

j  j  1=1  jl  j  jl  j 

where  AO,  A1,  C  and  D  are  the  model  coefficients,  is  the 
path  bearing  angle  measured  positive  clockwise  from  north  at 
the  jth  (  W,X,Y  or  Bj  station  and  T  is  the  path  range  to  the 
jth  station.  The  GRB  model  is  relatively  complex  and  is 
expected  to  exhibit  superior  performance.  The  extensive 
mcdel  is  based  on  knowing  the  distance  overland  (TL) ,  the 
distance  over  water  (TS) ,  the  total  distance  (T) ,  and  the 
path  bearing  angle  3  (Figure  2.5). 

It  was  noted  that  the  calibrated  mcdel  was  expected 
to  be  accurate  and  applicable  only  over  the  extent  of  ranges 
and  bearing  angles  embodied  in  the  calibration  data.  Hence, 
outside  the  region  covered  by  the  calibration  data  the  model 
may  not  be  as  accurate  as  within  the  data  coverage  region. 
Using  the  GRB  model,  a  drms  value  of  approximately  0.8  psec 
was  expected  in  areas  where  land  data  alone  was  used  to 
calibrate  the  model.  Inclusion  of  sea  calibration  data 
produced  a  drms  value  of  3.  35  to  0.50  jisec. 
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where 
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T  - -  psec, 
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psec. 


=  ith  secondary  stat ion-to-user  great-circle  path  length 


=  |heg^|ee^of  ag#al ectromagnetic  wave  in  a  vacuua 

=  surface  refractive  index 
=  1.000338, 

«  master  station-to- user  great-circle  path  length. 


=  coding  delay  found  in  Table  II, 

=  TD  bias  associated  with  the  ith  secondary 
station  (psec)  (Table  III)  . 

TABLE  III 


r 


TD  Bias  (b)  -  psec 


TD  W  -0.854 

TDX  -1.173 

TD  Y  -0.353 


SF  =  0.5  (-S1  ♦  S2  ♦  S3  -  S4  ♦  S5  ♦  S6)  (2.9) 

2. 

The  term  SI  is  the  SF  of  a  land  path  of  length  T 

s  j 

(psec)  from  the  jth  station: 

SI  =  0.795/T  ♦  0.439  ♦  (0.00245)  T  usee.  (2.10) 

sj  si 

The  terms  S2,  S3,  and  S4  combine  to  make  up  the 

secondary  factor  for  the  seawater  path  lengths. 
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S2  is  the  SP  for  the  total  path.  S3  is  the  SP 
using  the  seawater  coefficients  for  the  portion 
with  seawater,  and  S4  is  the  SP  using  the  seawater 
coefficients  for  the  land  path  distance.  The 
seawater  coefficients  are  found  in  Table  IV. 

TABLE  IV 

Seawater  Coefficients 


m 

128.8 

al 

3 

a2 

3 

0.  187 

a3 

* 

0.000652 

bl 

- 

3.188 

b2 

3 

-0.594 

b3 

3 

0.000329 

S?  (T)  =  b  1/T  ♦  b2  +  (b3)  T  p sec,  (2.11a) 

if  10  S  T  <  540  pssc, 

SF  (T)  =  a  1/T  ♦  a2  ♦  (a3)  T  ;isec,  (2.11b) 

if  T  >  540  psec. 

Term  S5  and  S6  are  the  SFs  of  land  paths  of  length 
T  and  T  .  SF5  is  the  Secondary  Factor  for  the 

j  lj 

total  length  using  the  land  coefficients  whereas 
SF6  is  the  Secondary  Factor  for  the  distance  over 
land  using  the  land  coefficients.  The  land 
coefficients  are  found  in  Table  V. 


TABLE  V 


Land  Coefficients 


AO 

M 

1.428 

D 

x2 

C 
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y2 
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m2 

-0.355 

2.  DMAHTC  MODEL 

1 .  ?ea  SF  Model 

The  equations  for  the  Sea  SF  Model  is: 


SF  =  (B1/T)  +  B2  ♦  (B3  T)  usee,  if  13  <  T  <  537  /isec,  (2.13a) 
SF  =  (A1/T)  ♦  A2  ♦  (A3  T)  usee,  if  T  >  537  /isec,  (2.13b) 

where  T  is  the  primary  phase  delay  (or  range)  in  microse¬ 
conds  (usee);  Ak  and  Bk  (k  =  1 ,  2,  and  3)  are  the  sea  model 
coefficients  used  by  DMAHTC  in  program  TDGSID  (Funakoshi, 
1982  ].  The  coefficients  are  found  in  Table  VI. 

2.  Land  SF  Mo^el 

The  solution  used  to  resolve  DMAHTC  Tables  ASF 
Corrections  is  called  Millington's  He-.hod  [DMAHTC,  1981]. 
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TABLE  VI 


Sea  SF  Model  Coefficients 


A  1 

3 

129.04323 

A  2 

M 

-0.40758 

A3 

S 

0.00064576813 

B 1 

—» 

2.741282 

B  2 

3 

-  0. 011402 

B3 

= 

0.00032774815 

This  method  is  based  on  the  premise  that  the  phase  distor- 
tion  due  to  a  composite  land-sea  path  is  the  arithmetic 
average  of  the  phase  distortion  found  in  the  forvard  and 
reverse  paths  of  the  propagated  signal  [  DMAHTC,  1981].  For 
example,  in  Figure  2.6  two  azimuths  have  been  drawn  on  the 
map  and  are  labeled  as  210°  and  23  5°.  Also  placed  on  the 
map  are  the  proper  ground  conductivities  which  have  been 
determined  in  the  field  by  the  Coast  Guard  calibration  team. 
A  great  circle  drawn  on  the  appropriate  chart  or  charts  from 
the  LORAN-C  Station  coordinates  to  the  area  under  considera- 
xion  spans  various  lengths  of  land  and  seawater.  Each 
length  or  segment  will  have  a  specific  conductivity  and 
distance.  The  total  of  these  conductivity  segments,  consti¬ 
tuting  a  land-water  profile,  will  determine  the  total  phase 
retardation  along  xhat  path  [U.S.  Naval  Oceanographic 
Office,  1  982  ].  All  azimuths  and  distances  are  computed 
based  on  the  World  Geodetic  System  (WGS)  datum  [DMAHTC, 
1981].  The  values  of  phase  retardation  for  a  given  ground 
conductivity  are  tabulated  in  the  National  Bureau  of 
Standards  (N3S)  Circular  57  3  [Speight,  1982]. 

The  formula  used  to  derive  the  ASF  Correction  for 
the  time  difference  for  a  master-slave  transmitting  station 
pal  r  is : 

ASF  Correction  =  (-Slave  Error)  -  (-Master  Error) 
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Figure  2.6  Composite  Land-Sea  Path 

The  mean  values  derived  for  one  station  from 
forward  and  reverse  solution  of  Millington's  Method 
subtracted  from  the  Sea  SF  Model.  The  differences 


presented  in  an  azimuth  array.  This  array  is  a  series  of 
geodetic  aziauths  radiating  froa  the  transaitter  with 
corrections  coaputed  at  incremented  distances  along  each 
azimuth.  After  the  values  are  computed  in  the  azimuthal 
array  a  computer  program  rearranges  them  into  a  aatrix  fora. 
The  aatriz  form  is  the  arrangement  of  corrections  into  rows 
and  columns  covering  a  specified  geographical  area  at  a 
constant  spacing.  The  purpose  of  the  matrix  is  to  enable 
the  corrections  froa  two  LORAN-C  transmitters  to  be  added 
algebraically ,  combined  into  a  single  matrix,  and  arranged 
in  the  desired  tabular  fora.  This  tabular  fora  is  the  body 
of  the  table  fDHAHTC,  1981]. 

3 .  J&b  le  Description 

Each  table  contains  a  complete  chain.  Figures  2.5, 
2.6,  and  2.7  depict  LORAN-C  ASP  Correctors  for  chains 
994  0-N,  X,  and  I  for  Honterey  Bay,  California.  A  table 
section  is  prepared  for  each  station  pair  (master  station 
and  one  slave  station)  in  a  LORAN-C  chain.  As  a  rule  the 
limits  of  the  table  coverage  are  determined  by  the  range  of 
the  ground wave  transmissions  for  the  LORAN-C  chain.  Each 
page  of  corrections  in  the  table  covers  an  area  three 
degrees  in  latitude  by  one  degree  cf  longitude,  with  correc¬ 
tions  printed  in  increments  of  five  minutes  of  arc.  Rate 
designation  and  page  numbers  are  printed  at  the  top  of  each 
correction  page.  Those  pages  where  latitude  and  longitude 
limits  contain  both  land  and  sea  are  included  but  correc¬ 
tions  apply  only  for  the  area  covered  by  the  U.S.  Coastal 
Confluence  Zone  (CCZ)  .  Large  land  bodies  and  areas  ouside 
the  CCZ  are  represented  by  blank  spaces  on  the  page.  ASF 
correction  values  can  be  either  positive  or  negative  (posi¬ 
tive  values  are  shown  without  sign) .  Areas  requiring  no 
correction  show  a  zero  value  which  in  some  cases  in  preceded 
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Figure  2.7  9940-W  ASF  Correction  Table 


by  a  negative  sign  indicating  that  the  zero  results  from  the 
rounding  off  of  a  value  slightly  less  than  zero  (indicates 
the  trend  of  the  correction) . 
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LONGITUDE  HEST 


Figure  2.8  9940-X  1SF  Correction  Taoie 


The  table  can  be  entered  directly  by  using  the 
ship's  position  determined  to  the  nearest  five  minutes  of 
arc  in  latitude  and  longitude  either  by  dead  reckoning  or 
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Figure  2.9  9940-1  ASF  Correction  Table 


some  other  means.  To  find  the  page  with  the  appropriate 
correction,  the  Page  Indexes  of  the  table  should  be 
utilized.  These  indexes  show  the  limits  and  page  number  of 
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all  pages  in  the  table.  To  locate  the  number  ◦£  the  page  on 
which  the  desired  correction  is  to  be  found  the  Page  Index 
is  entered  with  the  ship' s  position.  In  some  cases  the 
ship's  position  will  fall  on  the  page  limit  ir  either  lati¬ 
tude  or  longitude  or  both.  These  positions  are  repeated  on 
both  pages  and  either  page  aay  be  used. 

The  ASF  Correction  is  added  algebraically  to  the 
tiae  difference  for  the  LORAN-C  pair.  Interpolation  of  this 
data  will  not  necessarily  ia prove  the  accuracy  due  to  the 
method  used  to  determine  ASF  Correctors  [DMAHTC,  1981]. 
Since  the  correctors  are  computed  in  the  azimuthal  array  and 
are  based  on  the  conductivity  and  distance  over  which  the 
LORAN-C  electromagnetic  wave  travels,  the  ASF  Corrector 
between  the  published  ASF  Correctors  in  the  tables  may  not 
be  the  linear  interpolated  values.  For  example,  the  ASF 
Corrector  for  a  distance  of  500  m  with  an  azimuth  of  180°  is 
equal  to  1.5  psec.  The  ASF  corrector  for  a  distance  of  500 
m  with  an  azimuth  of  18 1  o  is  equal  to  1.6  jisec.  The  inter¬ 
polated  value  between  1.5  and  1.6  is  1.55.  The  true  ASF 
Corrector  for  the  180°  30'  azimuth  is  1.4  since  the  land 
distance  for  the  same  azimuth  is  less  than  the  land  distance 
for  the  180°  and  181®  azimuth.  The  LORAN-C  signal  passed 
over  a  harbor  [Dansford,  19  82  ]. 

F.  ATTEMPTED  DETERMINATION  OF  ASF  CORRECTORS  BY  HYDRO  FIELD 
PARTIES 

One  of  the  major  problems  encountered  by  hydrographic 
survey  operating  units  when  using  LORAN-C  for  position 
control  is  the  determination  of  the  ASF  Correctors  for  the 
survey  area.  The  four  surveys  mentioned  in  the  Introduction 
all  made  attempts  to  determine  the  correctors  by  comparing 
the  LORAN-C  rates  tc  a  second  source. 
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Bathymetric  Surveys  H-9822  [  NOAA  H-9822,  1979]  ana 

R-9823  [NOAA  H-9823,  1979  ]  Gulf  of  Alaska,  compared  the 

rates  from  an  Internav  LC-204  LORAN  receiver  to  computed 
rates  from  a  position  obtained  from  a  JHR-1  Satellite 
Navigation  Receiver  vhen  available.  Shore  ties  using  radar 
ranges,  visual  bearings,  and  sextant  angles  in  comparison  to 
LORAN-C  rates  were  also  made  prior  to  and  after  each  survey. 
The  calibrations  of  LORAN-C  rates  were  based  on  the  satel¬ 
lite  positions  only  since  the  positional  computation  of 
LORAN-C  and  JHR  Doppler  Satellite  were  made  on  the  WGS  1972 
datum  whereas  the  land  ties  were  based  on  the  NAD  1927 
datum. 

Bathymetric  survey  SO— 40-7-82  which  extended  along  the 
Washington,  Oregon,  and  California  Coasts  used  LOPAN-C  as 
navigational  control.  LORAN-C  time  differences  were 

compared  with  SATNAV  positions.  The  report  did  not  indicate 
whether  any  correctors  were  applied  [ NO AA  Ship  Surveyor. 
1982]. 

Finally,  Hydrographic  Survey  S-K902-Wh-82, 

Reconnaissance  Survey  of  Safety  Fairways,  Gulf  of  Mexico 
used  LORAN-C  as  a  positioning  control  system.  The  posi¬ 
tioning  unit  was  an  LC-204  receiver.  LORAN-C  rates  were 
input  via  the  HIDROPLOT  Controller,  a  special  purpose 
input-output  interface  which  is  the  nucleus  of  the  computer 
system  hardware  [Ombach,  1976  ].  Positions  were  computed  and 
plotted  by  Program  RK121 ,  LORAN-C  Real-Time  HYDROPLOT 
[  Backus ,  1  980  ]. 

ASF  Correctors  for  LORAN-C  were  achieved  by  visual  cali¬ 
bration  using  three  point  sextant  fixes  using  charted  oil 

« 

rigs  as  control  in  the  vicinity  of  the  survey  area.  A 
three-point  sextant  fix  is  a  convenient  and  accurate  method 
for  determining  the  position  of  a  hydrographic  survey 
vessel.  Sextants  are  used  to  measure  two  angles  between 
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three  objects  of  known  geographic  positon.  The  center 
object  is  common  to  both  angles.  The  position  of  the  obser¬ 
vers  taking  the  angles  is  fired  by  the  intersection  of  three 
circular  lines  of  position  [Ombach,  1976]. 

These  sextant  angles  were  recorded  and  later  tranferred 
onto  their  respective  charts  using  a  plastic  three-arm 
protractor.  A  plastic  three-arm  protractor  is  transparent 
and  made  up  of  one  fixed  arm  and  two  movable  arms  which 
contain  a r.  etched  line  that  is  radial  with  the  center  of  the 
protractor  [Ombach,  1976  ].  Sextant  angles  observed  in  pairs 
for  a  resection  fix  with  a  common  center  mark  may  be  plotted 
directly  by  this  instrument.  When  the  three  arms  are  placed 
at  the  angles  observed  and  fitted  so  as  to  pass  through  the 
plotted  positions  of  the  observed  stations  on  the  field 
sheet,  the  hole  at  the  center  of  the  three-ar^  protractor  is 
the  fixed  positon  of  the  vessel  [Ingham,  1975], 

Partial  correctors  for  each  area  surveyed  were  defined 
by  comparing  the  observed  rates  and  the  determined  rates 
plotted  on  the  nautical  chart.  The  partial  correctors  were 
applied  via  the  HYDRO PLOT  Controller.  However,  even  after 
applying  these  correctors,  the  plotted  position  still  disa¬ 
greed  with  the  ship's  determined  postion  with  respect  to  the 
oil  rigs.  Ship's  personnel  attributed  the  discrepancies  to 
one  or  more  of  the  following: 

1)  Accuracy  of  the  charted  rigs, 

2)  Weather  effect  on  LORAN-C, 

3)  Time  of  day, 

4)  Propagation  of  signal  over  land  path, 

5)  Three-arm  protractor  accuracy,  and 

6)  Error  in  the  conversion  by  the  software 
of  the  LORAN-C  rates  to  latitude  and 
longitude  [NOAA,  1982]. 
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There  is  an  apparent  need  for  a  LORAN-C  calibration 
routine  aboard  NOAA  ships  which  provides  the  ASP  Correctors 
for  program  RK121,  LORAN-C  Real-Tine  HYDRO  PLOT.  The  routine 
should  use  the  sane  geodetic  distance  computation  found  in 
RK121  and  use  the  sane  datum  as  that  of  the  naurical  chart 
of  the  survey  area. 

The  above  mentioned  discrepancies  illustrate  the  defi¬ 
ciencies  in  applying  a  single  ASF  Corrector  to  LORAN-C  data. 
The  accuracies  for  hydrography  cannot  be  met  using  single 
correctors  because  the  errors  are  non-linear  and  systematic. 
They  cannot  be  distributed  like  residuals  in  a  traverse. 
Schnebele  [  1979  ]  has  already  proven  that  single  ASF  area 
correctors  to  LORAN-C  positions  do  not  meet  the  accuracy 
standards  of  the  NOS  Hydrographic  Manual. 

Based  on  visual  inspection  of  the  DMAHTC  LORAN-C 
Correction  Tables,  ASF  Correctors  should  be  updated  every 
five  minutes  of  Latitude  or  Longitude  change.  In  Monterey 
Bay,  California,  there  is  approximately  0.1  to  0.2  jisec 
difference  for  every  five  minutes  of  change,  a  potential 
error  of  55  to  110  meters. 


4 


4 


V*  ‘m—  *r  J* — 


III.  EXPERIMENT II  PROCEDURE 

A.  FIELD  PROCEDURES 

In  order  to  compare  the  use  of  differential  LOHAN-C  with 
ASF  multiple  correctors,  typical  survey  operations  were 
planned  for  the  southern  portion  of  Monterey  Bay, 
California.  This  survey  was  conducted  in  conjunction  with  a 
comparative  evaluation  of  multiple  lines  of  positon  for 
selective  positioning  methods  [Anderson,  1982].  Four  micro- 
wave  ranging  systems  were  set  on  known  geographic  positions 

TABLE  VII 

Geographic  Names  and  Positions 


Microwave  System 
Stations 

Geographic  Position 
(NAD  1927) 

Dates 

1 

L  us&d. 

Seaside  4  (1964) 

360 

36* 

23.44596" 

June 

3-5,1982 

121® 

51' 

38.83281" 

Use  Mon  Ecc. 

360 

36* 

04.73031" 

June 

3-5,1982 

121  o 

52' 

35.98040" 

Geoceiver  Ecc. 

360 

36* 

32.49281" 

June 

3-5,  1982 

121  o 

53* 

25.21162" 

Mussel  Ecc. 

36  o 

36' 

18.25484" 

June 

3-5, 1982 

1210 

54 ' 

11.49661" 

Park  (1931) 

36  o 

53' 

13.90600" 

June 

6-7,1982 

121  o 

49' 

46.74300" 

Mulligan  RM1 

360 

44* 

56.49531" 

June 

6-7,1982 

121  o 

47* 

52.31090" 

Range  7  (1972) 

360 

39* 

02.47787" 

June 

6-7,1982 

1210 

49* 

08. 5  82  02" 

Mussel  (1  932) 

36  o 

37' 

18.15100" 

June 

6-7, 1982 

1210 

54' 

11.49661" 
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listed  in  Table  VII.  A  series  of  tracklines  were  ran  in  tvo 
separate  areas  as  shown  on  Figure  3.1.  To  ensure  that  the 
oicrovave  postioning  sys tea  was  working  properly,  the  equip¬ 
ment  was  calibrated  over  known  baselines  of  1497.47  aeters 
and  7877.31  aeters  at  the  beginning  and  end  of  the  project. 
Trackline  observations  were  only  made  during  the  daytime  in 
fair  weather  conditions  so  as  to  eliminate  sky  waves  and 
weather  changes  that  influence  LORAN-C  signal  propagation 
characteristics  [Samaddar,  1982  ].  The  vessel  used  was  the 
126  foot  R/V  Acania  which  is  operated  by  the  Naval 
Postgraduate  School. 

The  positioning  equipment  consisted  of  a  Micrologic 
ML- 1000  LORAN-C  receiver  (0.01  usee  resolution)  and  a 
Trisponder  Microwave  System  provided  by  Racal-DECCA  Survey, 
Inc.  The  Trisponder  Microwave  System  consisted  of  four  DNT1 
Model  21 7c  transponders,  four  DNT1  Model  21017  HP  sector 
antennas  with  87°  by  5°  beam  widths,  one  DNT1  DDMU  (Digital 
Distance  Measuring  Unit),  two  Omni  DVTI  Model  21019  HP 
antennas,  a  Houston  Instruments  Model  DP3-M2D/RC3  plotter, 
and  a  Texas  Instruments  743  terminal  (Table  VIII).  The 
manufacturer's  published  accuracy  for  the  positioning  equip¬ 
ment  is  t  1  m  for  a  single  range  [Racal-DECCA  Survey,  1981]. 
Anderson  [  1982  ]  discusses  the  accuracy  of  four  lines  of 
position.  The  four  Decca  Trisponder  distances  were  recorded 
via  a  Texas  Instruments  743  data  terminal  while  the  LORAN-C 
rates  were  manually  logged.  The  data  was  acquired  at  one 
minute  intervals  while  the  ship  maintained  constant  course 
and  speed.  The  recorded  LORAN-C  rates  were  9940-W,  X,  and  Y 
of  the  West  Coast  chain. 

To  test  the  potential  for  calibrating  the  LORAN-C  System 
using  the  Semi- Empirical  Model,  the  correction  tables,  and 
multiple  observed  field  correctors,  the  positions  derived 
from  the  microwave  system  measurements  were  used  to  compute 
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expected  LORAN-C  tine  differences  at  each  point.  The 
difference  or  offset  between  these  expected  time  differences 
and  the  observed  values  were  computed  for  the  three  methods. 
The  mean  offset,  standard  error,  and  drms  values  were  also 
computed  and  compared. 

B.  MICROWAVE  SYSTEM  POSITIONING 

The  geographic  position  of  the  ship  based  on  four  lines 
of  position  was  determined  using  a  computer  program  called 
GP3YLQ  (Geographic  Position  by  Least  Squares)  written  by  the 
author  (See  Appendix  A)  .  GPBYLQ  contains  subroutine  LSQR 
(Least  Squares),  which  Is  a  least  squares  adjustment  written 
by  Paul  R.  Wolf,  Ph.D.  [1974]  and  revised  by  LCDR  D.  Leath 
[1981].  Geographic  postions  were  converted  to  X,Y  (meters) 
which  in  turn  were  converted  to  geographic  position  via 
subroutines  GPTOXY  and  XYTOGP,  respectively  [Wallace,  1974  ]. 


Subroutine  GPTOXT  and  XTTOGP  are  based  on  the  Bodified 
Transverse  Mercator  Grid  which  was  centered  in  the  survey 
area. 

The  Modified  Transverse  Mercator  (MTM)  projection  is 
used  by  the  National  Ocean  Survey  and  is  sinilar  to  the 
projection  used  in  the  Universal  Transverse  Mercator  (UTB) 
systea.  The  aain  difference  is  that  in  the  MTM  a  Central 
Meridian  is  picked  that  is  near  the  survey  area  instead  of 
being  fixed  at  a  particular  meridian  [Wallace,  1971]. 
Central  Meridian  (CMER)  ,  False  Easting  (FEST)  ,  and 
Controlling  Latitude  (CLAT)  are  the  three  paraaeters  which 
define  the  MTM  projection.  CMER  is  the  aean  longitude 
coaputed  using  the  maxiaua  and  ainiaua  longitudes  of  the 
survey  limits,  FEST  is  the  X-Coordinate  that  is  assigned  to 
the  Central  Meridian,  and  CLAT  is  the  distance  in  aeters 
from  the  equator  to  seme  reference  latitude  [Wallace,  1971]. 
The  Central  Meridian,  False  Easting,  and  Controlling 
Latitude  used  for  Monterey  Bay,  California  referenced  to  SOS 
Chart  18685  are: 

CMER  =  1210  56'  00.0", 

FEST  =  20000.0, 

CLAT  =  4050000.0, 

To  be  consistent  with  the  National  O'-ean  Survey  charts 
of  the  area,  all  computations  were  done  relative  to  North 
American  Datum  (NAD)  19  27  geographic  positions.  All 
programs  were  executed  on  an  IBM  30  33  computer  located  at 
W.R.  Church  Computer  Center,  Naval  Postgraduate  School, 
Monterey,  California. 

C.  LORAN-C  COMPUTATIONS 

The  differences  or  offset  between  the  observed  and 
computed  LORAN-C  rates  using  the  Semi-Empirical  TD  Model, 
ASF  LORAN-C  Correction  Tables,  or  the  Multiple  Observed 
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Correctors  were  compared  to  the  offsets  between  the  observed 
TD  rates  and  the  computed  rates  for  which  only  the  seawater 
Secondary  Factors  (SF)  were  applied.  The  comparison  of  the 
offsets  between  the  four  methods  illustrates  the  improvement 
in  positional  accuracy  after  applying  ASF  Correctors. 

1  •  Seawater  Secondary  Factors  (SF) 

Tim's  differences  using  only  the  seawater  Secondary 
Factors  for  each  of  the  geographic  positions  were  computed 
using  program  L08AN  written  by  the  author  (Appendix  C) . 
Seawater  Secondary  Factors  (SF)  were  computed  using  formula 
2.13  and  the  coefficents  found  in  Table  71.  All  TOA 
distances  in  meters  were  determined  using  subroutine  INVER1. 
INVER1  is  a  geodetic  inverse  routine  using  I.  vincenty's 
modified  Rainsford's  method  with  Helmert's  elliptical  terms, 
programed  by  LCDR  L.  Pfeifer,  NOA  A  [  1975  ].  Subroutine 
INVER  1  is  accurate  to  0.0  001  m  halfway  around  the  world 
[Pfeifer,  1982  ].  All  distances  were  converted  to  microse¬ 
conds  using  299.792458  m/fisec.  Time  differences  (TD)  were 


computed  fro*  equation  2.5.  North  American  Datum  1927 
geographic  positions  were  used  for  all  computations.  See 
Table  IX  for  the  positions  of  LORAN-C  9940  transmitters 
[Hiordan,  1979]. 


To  ensure  that  subroutine  INVER1  was  functioning 
properly  the  distances  between  the  master  and  secondary 
stations  were  compared  to  the  NOS  published  baseline 
distance  [Riordan,  1979  ].  The  published  distances  and  the 
results  from  routine  INVER1  are  listed  below. 

Published  Baseline  Computed  Baseline 

Distance  (m)  Distance  -  INVERT  (m) 


9940-H  837,777.  0929 
9940-X  327,886.  3720 
9940-Y  589,298.  5712 


837,777. 115 
327,886.  316 
589,298.589 


The  difference  between  the  published  and  computed  baselines 
ranged  from  0.02  to  0.06  m. 

Differences  or  offsets  (x t  )  were  obtainad  by 
subtracting  the  observed  LORAN-C  rates  from  the  computed 
values  from  the  various  methods.  The  mean  difference  or 
offset  (x)  and  standard  deviation  (s)  in  microseconds  for 
each  rate  were  determined  using  equations  3.1  and  3.2 
[Wonnacott,  1935  ]: 


Mean  offset  (/isec) : 
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where:  x  =  original  observation  in  usee, 
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The  mean  offset  in  microseconds  can  be  converted  to  aeters 
using  equation  3.3  [Heinzen,  1977]: 


Mean  offset 
(aeters)  : 

where  w 


z  v 


sin  a 


(3.3) 


distance  corresponding  to  one  aicrosecond 
on  the  baseline 

149.896229  eaters  [Bigelow,  1963]. 

onerhalf  the  angle  between  the 
radius  vectors  from  the  position  to 
the  master  and  secondary  stations. 


Equation  3.4  was  used  to  coapute  drms  values  in  aeters 
[Heinzen,  1977]: 
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where:  p  =  correlation  coefficient  =  0.33, 

s  and  a  are  a?  above  with  the  subscript 
denoting  the  appropriate 
secondary  station. 

The  correlation  coefficient  (p)  is  a  result  of  the  secondary 
station  having  a  common  line  of  position  with  the  master 
station.  Although  often  ignored,  various  authors  assign 
values  ranging  from  0.33  to  0.40.  Bigelow  [  1963]  chooses 
P  =  0.33. 

2.  Semi- Empirical  TD  Srid 

To  determine  if  the  Semi-Empirical  TD  Grid  would 
reduce  either  the  drms  value  or  the  offset:  between  the 
observed  and  expected  TD  rates,  program  LOPLC  (Line  of 
Position  -  LORAN-C)  was  written  by  the  author  (Appendix  B)  . 
Program  LOPLC  computes  the  distance  over  land,  distance  over 
water,  the  total  distance,  and  the  azimuth  from  north  for 
the  transmitting  station  using  Subroutine  INVER1.  The 
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land-sea  distances  were  cos pu ted  by  selecting  a  point  which 
was  located  along  the  coast  of  Monterey  Bay.  Subroutine 
SHOBPT  (Shore  Point)  interpolated  &  geodetic  position  from 
23  geodetic  points  which  outlined  Monterey  Bay,  California. 
The  land/sea  Secondary  Factor  was  detemined  in  Subroutine 
SECFAC  which  is  found  in  progras  LOPLC  using  equations  2.8 
through  2.12.  All  computations  ware  based  on  MAO  1927 
geodetic  datum. 

Tests  were  made  on  program  LOPLC  using  data  found  in 
the  Semi- Empirical  TO  Grid  article  [The  Analytic  Science 
Corporation,  1979]  using  the  MGS  1972  datum.  Station  TASC 
55  located  at  latitude  34°  34'  18.3"  E  and  longitude  120° 

39'  40.3"  W,  was  selected  from  The  Analytic  Science 

Corporation  article.  It  was  one  of  the  stations  used  to 
calibrate  the  coefficients  for  the  Semi-Empirical  TO  equa¬ 
tions  discussed  earlier.  The  only  distances  listed  for  TASC 
55  were  the  individual  distances  over  land  and  over  water. 
The  total  distance  between  the  transmitters  and  TASC  55  was 
computed  by  adding  the  land  and  sea  distances.  The 
following  station-to-sits  path  segment  lengths  are  listed 
for  TASC  55: 


Land 

Distance 

Sea 

Dist  ance 

TASC  55 

Total 

Dist  ance 

IMVER1 

Total 

Distance 

Station 

(111) 

(Isa) 

(Isa) 

(Isa) 

Master 

540.730 

35.548 

576. 278 

576.08  3 

X-Secondary 

369.248 

126 .096 

495.344 

495. 178 

Y-Seccndary 

525.659 

15 .832 

541.491 

541.308 

Using  the  land  distances,  sea  distance,  and  total  distances 
to  TASC  55,  and  the  computed  atimth  from  subroutine  INVEP1 , 
program  LOPLC  produced  offsets  of  3.86  ;isec  and  1.13  /isec 
for  9940-X  and  9940-1  respectively  when  compared  to  the 
expected  time  differences  at  TASC  55. 


An  attempt  was  made  to  determine  if  this  discrepency 
was  due  to  program  LOPLC.  &  comparison  of  the  total 
distances  in  the  report  between  the  transmitters  and  TASC  55 
were  compared  to  the  total  distances  computed  by  INVEE1. 
The  difference  in  distance  between  Subroutine  INVER1  and  the 
total  distances  from  the  Semi-Empirical  TD  Srid  acticle 
ranged  from  160  -  200  m  corresponding  to  0.5  to  0.65  usee. 
This  results  in  time  difference  errors  of  0.04  usee  for  the 
9940-X  rate  and  a  0.  10  jusec  for  the  9940-Y  rate,  implying 
that  the  offsets,  0.86  jasec  and  1.13  jisec,  are  caused  in 
part  by  the  method  in  which  the  azimuth  from  north  is  deter¬ 
mined.  Unfortunately,  azimuth  data  from  TASC  55  was  not 
presented  in  the  article.  The  Analytic  Science  Corporation 
has  been  contacted  on  numerous  occasions  in  an  attempt  to 
ascertain  their  method  of  determining  distance  and  azimuth. 
As  of  this  date  there  has  been  no  response. 

Nevertheless,  the  data  from  Monterey  Bay  was 
utilized  in  Program  LOPLC  to  obtain  results  that  could  be 
compared  to  that  obtained  by  the  other  methods.  If  this 
method  is  accurate  enough,  ASF  Correctors  could  be  deter¬ 
mined  via  computer  for  each  individual  position  without 
using  tables  or  field  determined  correctors.  The  mean 
offset  (x)  and  standard  error  (s)  in  microseconds  between 
the  observed  and  calculated  rates  were  computed  using  equa¬ 
tions  3.1  and  3.2,  respectively.  The  mean  offset  (1)  and 
drms  in  meters  were  computed  using  equation  3.3  and  3.4, 
respectively. 

3.  Calculated  Ta ble  ASF  Correctors 

The  offset  between  the  observed  LORAN-C  rates  and 
the  expected  values  with  applied  ASF  Correctors  from  the 
LORAN-C  Correction  Tables  and  the  seawater  Secondary  Factor 
(equation  2.13)  were  determined  using  program  LORTAB  which 


was  written  by  the  author  (Appendix  D)  .  The  ASF  Corrector 
for  each  position  was  selected  using  subroutine  TABLE. 
Subroutine  TABLE,  which  is  found  in  program  LORTAB  deter¬ 
mines  an  ASF  Corrector  for  each  data  point  based  on  its 
geodetic  position.  The  ASF  Correctors  used  in  subroutine 
TABLE  (see  Figures  2.5,  2.6,  and  2.7)  are  located  between 

latitude  36°  35'  N  and  36<>  55'N  and  longitude  1210  50'  W  and 
122°  00'  W.  The  difference  in  sign  between  the  ASF 

Correctors  in  subroutine  TABLE  and  those  found  in  the 
LORAN-C  Correction  Table  is  due  to  the  difference  in  their 
application.  ASF  Correctors  from  the  tables  are  applied  to 
observed  rates  while  ASF  Correctors  from  Subroutine  TABLE 
are  applied  to  the  calculated  LORAN-C  rates.  Negative  ASF 
Correctors  from  the  LORAN-C  Correction  Table  were  applied  to 
the  calculated  time  differences  to  be  consistent  with  the 
application  of  Secondary  Factors  to  the  computed  primary 
phase  delay. 

As  before  equations  3.1  through  3.4  were  used  to 
compute  the  the  mean  offset ,  standard  error,  and  drms.  If 
this  application  of  LORAN-C  Correction  Tables  is  accurate 
enough,  it  precludes  the  need  to  determine  ASF  Correctors  in 
the  field. 

4 .  Observed  ASF  Correctors 

Observed  ASF  Correctors  were  determined  using 
Program  ASFSEL  (ASF  Selection)  which  was  written  by  the 
author  (Appendix  E)  ,  Schnebele's  prior  data  [  1979  ],  and  the 
June  1982  data.  Program  ASFSEL  (ASF  Selection)  was  written 
by  the  author.  This  program  computes  the  ASF  Correctors  by 
subtracting  the  observed  LORAN-C  rates  from  the  expected 
values.  Only  the  seawater  Secondary  Factors  from  equation 
2.13  have  been  applied  to  the  calculated  time  differences. 
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The  aean  ASF  Correctors  for  the  LORAN-C  rates  were  deter¬ 
mined  at  every  minute  of  latitude  and  longitude  between  36» 
50'  N  and  36°  35'  N  and  122°  04*»  and  122«  49*  W.  See 
Appendix  E  for  mean  ASF  Correctors  at  one  minute  intervals. 
The  ASF  Correctors  were  then  selected  and  assigned  to 
subroutine  TABLE  in  Program  LORTAB  at  five  minute  latitude 
and  longitude  intervals.  All  values  were  entered  to  the 
nearest  hundredth  of  a  microsecond.  Equations  3.1  through 
3.4  were  used  as  before  for  computations.  This  determina¬ 
tion  of  ASF  Correctors  in  the  field,  if  accurate  enough,  may 
allow  the  use  of  LORAN-C  as  a  positioning  system  for  hydro- 
graphic  surveying  in  the  future. 


IV.  RESULTS 


A  total  of  620  time  differences  (TD)  and  geographic 
positions  based  on  four  lines  of  position  were  recorded  in 
the  southern  portion  of  Montery  Bay  in  order  to  compare  the 
use  of  differential  LOBAN-C  with  ASF  Multiple  Correctors. 
The  data  was  divided  into  four  sets.  The  first  data  set  is 
Schnebele's  [1979]  prior  data  consisting  of  a  total  of  130 
data  points  collected  on  two  separate  days,  June  12  and  July 
25,  1979  between  latitude  36®  38*  N  and  36®  47'  N  and 
longitude  121®  49*  W  and  12  2°  02'  W.  The  recorded  LORAN-C 
time  difference  rates  were  9940-Y  and  9940-H. 

The  second  data  set  consists  of  193  time  differences 
collected  on  June  7,  1982.  The  recorded  LORAN-C  rates  were 
also  9940-Y  and  994 0-W.  Data  set  Number  2  is  located 
between  latitude  36°  40'  N  and  36®  45'  N  and  longitude  121® 
54'  W  and  122®  00'  W.  See  Figure  3.1  for  the  location  of 
test  areas  in  Monterey  Bay,  California. 

Data  set  Number  3  is  located  in  the  same  area  as  data 
set  Number  2,  between  latitude  36®  40'  N  and  36®  45'  N  and 
longitude  121®  W  54*  and  122®  00'  W.  This  set,  consisting 
of  128  points  with  recorded  rates  9940-X  and  9940-Y,  was 
recorded  on  June  6,  1982. 

The  final  set,  data  set  Number  4  was  recorded  between 
June  3  and  June  5,  1982.  It  contains  169  points  located 
between  latitude  36®  36  N  and  36®  39*  N  and  longitude  121® 
53*  W  and  121®  58*  W.  The  recorded  rates  were  9940  -X  and 
9940-Y  which  are  the  same  as  data  set  three.  The  data  for 
rates  9940-X  and  9940-Y  was  kept  separate  so  as  to  determine 
if  there  was  a  significant  difference  between  the  offshore 
(data  set  three)  and  inshore  (data  set  four)  drms  values  due 
to  phase  recovery  (see  Chapter  Two)  . 
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Number  of 
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L&t 
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Lon 

1 

June  12,  1979 
July  25,  1979 

9  940-1 
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122/02 
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June  7,  1982 

9  940-  Y 
9  940-W 
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36/40 

36/45 

121/54 

122/00 

3 

June  6,  198  2 

9  940- X 
9  940- Y 

128 

36/40 

36/45 

121/54 

122/00 

4 

June  5,  1982 

9  940-  X 
9  940- Y 

169 

36/36 

36/39 

121/53 

121/58 

Table  X  provides  a  convenient  breakdown  of  the  parame¬ 
ters  for  each  of  the  data  sets.  The  table  consists  of  the 
data  set  number,  the  date  the  data  was  collected,  the 
LOR AN -C  time  difference  rates,  the  number  of  data  points, 
and  the  area  limits  in  latitude  and  longitude.  It  defines 
the  parameters  for  the  four  data  sets  of  Tables  XI  (Seawater 
Secondary  Factor  Errors)  ,  XII  (Semi-Empirical  TD  Grid 
Correction  Errors),  XIII  (Multiple  LORAN-C  Correction  Table 
Errors)  ,  and  XIV  (Multiple  Observed  Correction  Errors)  .  All 
basic  data  and  computations  are  kept  on  file  with  the  NOAA 
hydrography  instructor  at  the  Naval  Postgraduate  School, 
Monterey,  California.  (Individuals  seeking  this  information 
should  contact  the  Oceanography  Department.) 

A.  SEAWATER  SECONDARY  FACTORS 

LORAN-C  rates  were  computed  using  only  the  seawater 
Secondary  Factor  (equation  2.13).  Ihe  LORAN-C  computed  time 
differences  are  basically  uncorrected  rates  since  no 
Additional  Secondary  Factors  were  applied.  Offsets  between 


TABLE  XI 

Seawater  Secondary  Factor  Errors 


Mean 

Mean 

Standard 

Offset 

Offset 

Error 
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Hates 
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-0.526 
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0.059 
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d) 
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X 

4 

(June  3-5,  1982,  169  data  points). 

1.582  6  50.089  r  0.055 

68.  1 

Y 

-0.817 

-237.109 

0.187 

the  observed  and  calculated  rates  using  only  the  seawater 
Secondary  Factor  were  computed  to  illustrate  the  improvement 
in  position  after  applying  the  Semi-Empirical  TD  Model,  ASF 
Loran-C  Correction  Tables,  or  the  Multiple  Observed  ASF 
Correctors.  The  mean  offset  in  microseconds  and  meters, 
standard  error,  and  the  drms  are  found  in  Table  XI. 
Examples  of  offsets  for  several  data  points  are  listed  in 
Appendix  B  after  program  LOBAN. 

It  was  stated  earlier  in  Chapter  One  that  Schnebele 
obtained  a  66  m  1  drms  using  Y  and  W  rates.  The  drms  of  66 
m  was  obtained  using  48  data  points  which  were  located  10  km 
or  more  offshore.  The  110.4  m  1  drms  for  the  Y  and  W  rates 
in  Table  XI  is  a  result  of  combining  the  130  inshore  and 
offshore  positions.  The  increase  from  66.0  m  1  drms  for  the 
offshore  positions  to  110.4  m  for  combined  offshore  and 
inshore  oositions  indicates  that  the  application  of  a  single 


ASF  Corrector  is  dependant  on  the  size  of  the  area.  The  48 
data  points  were  located  between  latitude  36°  41*  N  and  36° 
46*  N  and  longitude  121°  55'  W  and  122°  02*  W. 

Schnebele  also  obtained  the  66  i  1  dras  by  adjusting  the 
LORAN-C  observed  tiae  differences  which  were  shewed  due  to 
the  ship's  notion  and  the  five  second  averaging  interval  of 
the  LORAN-C  receiver.  These  caused  the  observed  TDs  to  be 
several  seconds  old  in  conparison  to  the  nicrowave  systea 
aeasurenents  [Schnebele,  1979].  Due  to  the  large  anounts  of 
data  fron  the  June  1982  survey  operations,  no  deskewing  was 
done. 

The  offset,  standard  error,  and  dras  for  data  set  Nunber 
2  (June  7,  1982)  compares  well  to  data  set  Number  1 
(Schnebele's  data  -  all) .  Also,  the  38.6  a  1  dras  for  the 
X-Y  rates  is  between  Schnebele's  42.0  a  1  dras  prediction 
[Schnebele,  1979]  and  Nelson's  findings  of  30  a  1  drms  in 
San  Francisco  Bay  [General  Electric  Co.,  1979].  The  large 
dras  of  68.1  a  for  data  set  Number  4  (X-Y  rates)  is  probably 
due  to  the  phase  recovery  of  the  electromagnetic  wave  from 
the  9940- Y  transmitter  located  in  Searchlight,  Nevada. 
Finally,  if  ASF  Correctors  are  not  applied,  dras  values 
ranged  from  38.6  a  for  9940  X-Y  to  101.4  m  for  9940  Y-W  for 
data  sets  2  and  3.  The  large  offsets  for  all  data  sets 
indicate  a  systematic  error,  ranging  from  150  m  to  700  a, 
which  precludes  the  use  of  this  method  for  hydrographic 
surveying. 

B.  SEHI -EMPIRICAL  TD  GRID 

To  determine  if  the  Semi-Empirical  TD  Grid  would  reduce  the 
offset  and  dras  for  hydrographic  surveying,  program  LOPLC 
(Appendix  C)  was  applied  to  data  sets  1  through  4.  The 
Semi- Empirical  TD  Grid  applies  a  Secondary  Factor  and  an 
Additional  Secondary  Factor  to  the  primary  phase  delay  based 


or  th«  distance  overland,  the  distance  over  water,  and  the 
total  distance  using  aean  land  and  water  conductivities. 
The  aean  offset  in  aeters  and  aicroseconds,  standard  error, 
and  dras  are  listed  in  Table  XIX.  Examples  of  offsets  for 
several  data  points  can  be  found  in  Appendix  C  after  prograa 

TABLE  XII 

Semi-Empirical  TD  Grid  Correction  Errors 


Mean 

Offset 

(lessee) 


Mean 

Offset 

(ill) 


Standard 

Error 

(s^isec) 


a)  Set  1  (Schnebele's  data,  130  data  points) 
Y  1.131  330.672  0.086 


1.145 


630.278 


0.  139 


b)  Set  2  (June  7,  1982  data,  193  data  points) 
Y  1.108  324.  975  0.077 

1.105 


"-a.  m  \  7  -J 

324.  975 
603.926 


gomti 
.077 
0.  126 


c)  Set  3  (June  6,  1982  data,  128  data 
X  1.131  386.  940 

0.589 


,  ua #  » 

386. 940 
172.923 


0.061 

0.084 


d)  Set  4  (June  3-5,  1982  data,  169  data  points) 


'1.056 

0.341 


uo.  ia  m  it. 

365.  1  13 
99.084 


0. 052 
0.  189 


dras 

(S) 


113.5 


112.7 


39.4 


68.2 


LOPLC.  Earlier,  Program  LOPLC  had  been  tested  with  data 
found  in  the  article  by  The  Analytic  Science  Corporation 
[  1979  ].  Program  LOPLC  produced  offsets  of  0.86  /isec  and 
1.13  psec  for  the  time  differences  from  rates  9940-X  and 
9940- Y,  respectively.  As  previously  stated,  the  large 
offsets  may  be  due  tc  the  method  by  which  the  distance  and 
azimuth  were  computed.  The  drms,  for  the  data  from  Monterey 
Bay,  obtained  with  the  Semi -Empirical  TD  Grid  was  similar  to 
the  drms  errors  for  seawater  Secondary  Factor  Model.  This 
might  be  an  indication  that  Program  LOPLC  is  correct  but  the 


bias  needs  to  be  adjusted  to  reduce  the  large  offset  between 
the  observed  and  calculated  TD  rates.  The  aean  offset 
ranged  from  99  to  630  ■.  Bias  reduction  could  be  achieved  by 
applying  land-sea  data  for  the  entire  West  Coast  to  the 
Se*i- Empirical  aodel.  Again,  the  existence  of  large  offset 
values  precludes  the  use  of  this  method  for  hydrographic 
surveying. 

C.  TABLE  ASP  CORHECTOBS 

To  deteraine  if  the  multiple  ASF  corrections  from  the 
LORAN-C  Correction  Table  would  reduce  the  offset  and  dras  to 
meet  the  NOS  accuracy  standards,  program  LORTAB  (Appendix  D) 
was  applied  to  data  sets  1  through  4.  ASF  Correctors  from 

TABLE  XIII 

Multiple  LORAN-C  Correction  Table  Errors 


Sates 


Mean 
Offset 
(x-pse  c) 


Mean 

Offset 

(III) 


Standard 

Error 

(s-^sec) 


drms 

(1) 


j  a)  Set  1  (Schnebele's  data,  130  data  points) 
Y  -0.229  -67.068  0.  104 

W  0.169  92.540  0.160 


b)  Set  2  (June  7.  1982  data.  192  data 


f  m  * 

-0.309 

0.123 


UCA  f  1  J  Cm 

-90.614 
S’7.  189 


Soints) 
.085 


0.085 
0.  144 


c)  Set  3  (June  6,  1982  data.  128  data 


I  ”  J  £  1 

0.§59 

-0.319 


ua  m  I  Cm\J 

191. 134 
-93.486 


Boints) 
.073 
0.  101 


d)  Set  4  (June  3-5,  1982  data,  169  data  points) 
Nc  ASF  Corrections  listed  in  Table  for  the 


123.6 


116.9 


47.3 


south¬ 


east  end  of  Monterey  Bay  next  to  the  shore  line. 


the  LORAN-C  Correction  Tables  are  determined  from  field 


observation  of  land  conductivities  by  the  U.S.  Coasn  Guard 


Calibration  Team  [Marina  Science  Department,  1982].  The 
results  are  listed  in  Table  XIII. 

When  coapared  to  the  Seavater  Secondary  Factor  Error  in 
Table  XI,  the  dras  values  using  LORAN-C  Correction  Tables 
were  increased  slightly  while  the  offsets  were  reduced 
substantially.  The  offsets  ranged  froa  65  -  200  a.  Since 
the  dras  ranged  froa  47.3  to  116.9  m  for  rates  9940  X-Y  and 
9940  Y-w  respectively,  the  application  of  ASF  Corrector  froa 
the  tables  does  not  meet  the  NOS  accuracy  standard.  Again, 
large  offsets  and  the  increase  in  dras  precludes  the  use  of 
this  aethod  for  hydrographic  surveying. 

D.  MOLTIPLE  OBSERVED  CORRECTORS 

To  deteraine  if  multiple  observed  correctors  would 
diminish  the  offset  and  drms  values,  mean  ASF  Correctors 
were  selected  at  one  minute  latitude  and  longitude  intervals 
using  Program  ASFSEL.  One  minute  ASF  Correctors  are  shown 
at  the  end  of  Program  ASFSEL  in  Appendix  E.  From  the  one 
minute  grid,  mean  ASF  Correctors  were  selected  and  entered 
into  Subroutine  TABLE  at  five  minute  latitude  and  longitude 
intervals  in  Program  LORTA3.  The  following  is  an  example  of 
the  9940-X  ASF  Correctors  at  five  minute  latitude  and  longi¬ 
tude  intervals  for  the  program: 


122/00/00.0 

121/55/00.0 

121/50/00.0 

36/50/00.0 

-1.  52 

-1.63 

36/45/00.0 

-1.  52 

-1.56 

-1.61 

36/40/00.0 

-1.58 

-1.60 

This  is  the 

same  format 

used  in  the 

LORAN-C  Ccrrec 

Tables.  Program  LORTAB  was  applied  to  the  four  data  sets. 
The  error  results  are  listed  in  Table  XIV. 

The  drms  value  obtained  with  multiple  observed  correc¬ 
tors  were  all  reduced  when  compared  to  the  drms  for  the 


TABLE  XIV 


Multiple  Observed  Correction  Errors 


Rates 


Mean  Mean  Standard 

Offset  Offset  Error 

(X-jisec)  !-■)  (s-^isec) 


drms 

(a) 


a) 

Set 

Y 

1 

(Schnebele's  data,  130 
0.033  9.585 

data  point: 
0.076 

W 

0. 

052 

2  8.604 

0. 

113 

b) 

Set 

2 

(June 

7.  1982, 

193  data 

points 

1)73 

Y 

0. 

0^8 

8.318 

0. 

w 

-0. 

016 

-8.944 

0. 

116 

c) 

Set 

3 

(June 

6.  1982, 

1 28  data 

points 

) 

X 

-0. 

0$4 

-11.490 

0. 

055 

Y 

0. 

010 

2.9  13 

0. 

066 

d) 

Set 

X 

4 

(June 

0. 

3.5,  1982 
001 

,  169  data  points) 
0.383  0.052 

Y 

-0. 

008 

-2.225 

0. 

187 

87.6 

89.3 

38.3 

67.5 


seawater  Secondary  Corrector.  The  most  impressive  reduction 
in  drms  was  within  Schnebele's  data  which  covered  an  area  of 
seven  minutes  of  latitude  and  12  minutes  of  longitude.  The 
drms  for  seawater  Secondary  Correctors  was  110.4  m  whereas 
the  drms  for  the  same  data  using  multiple  observed  correc¬ 
tors  was  87.6  m.  This  is  a  smaller  drms  than  that  of  the 
June  7,  1982  data  (data  set  Number  2)  which  was  obtained 
three  years  later.  It  appears  that  87.6  m  1  drms  is  nearly 
the  minimum  error  that  can  be  obtained  for  the  9940  Y-W 
rates  in  Monterey  Bay  after  applying  multiple  observed  ASF 
Correctors  at  five  minute  latitude  and  longitude  intervals. 
For  the  LORAN-C  rates  9940  X-Y,  a  38.3  m  1  drms  was  obtained 
for  same  five  minute  area  covered  by  data  set  Number  2  (June 
7,  1982).  LORAN-C  rates  9  940  X-Y  were  not  obtained  for  the 
same  size  area  covered  by  data  set  Number  1  (seven  minutes 
of  latitude  and  twelve  minutes  of  longitude  -  9940  Y-W)  due 
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to  the  restriction  of  ship  time  and  the  length  of  time  the 
Racal-DECCA  Trisponder  electronic  equipment  had  been  loaned. 

The  drns  for  the  June  3-5,  1982  inshore  data  was  only 
reduced  to  67.5  a  froa  68.1  a  for  seawater  Secondary 
Correctors.  The  small  change  in  error  at  the  coast  is 
probably  a  result  of  the  erratic  behavior  of  phase  recovery 
discussed  earlier  in  Chapter  Two.  The  drns  value  could 
presumably  be  reduced  if  the  correctors  were  applied  at  one 
ainute  intervals.  This  would  be  a  very  costly  method  of 
calibrating  Loran-C  for  hydrographic  surveying. 


7.  COMCIOSIOSS 


It  was  noted  in  Chapter  One  that  the  saallest  scale 
routinely  used  for  coastal  surveys  is  1  ;80,000.  This  yields 
an  allowable  error  of  40  m  1  drms  with  no  systenatic  errors. 
This  paper  deterained  whether  or  not  applying  aultiple 
Additional  Secondary  Factors  (ASF)  Correctors  to  LOBAN-C 
lines  of  position  would  reduce  the  dras  sufficiently  to  meet 
the  accuracy  standards  set  by  the  National  Ocean  Survey. 

Three  aethods  of  applying  aultiple  ASF  Correctors  were 
tested.  The  first  approach  coaputes  the  time  difference 
based  on  a  Semi-Eapirical  ID  Grid.  The  Semi-Empirical  Model 
produced  large  offsets  in  the  9940-B,  9940-X,  and  9940-Y 

time  differences.  The  offsets  ranged  from  99  to  630  m.  The 
drms  for  9940  X-Y  combination  was  39.4  m  and  the  drms  for 
9940  W-Y  combination  was  102.7  a. 

The  second  method  applies  AS?  Correctors  found  in  the 
DMAHTC  LORAN-C  Correction  Tables  to  LORAN-C  lines  of 
position.  The  application  of  the  tables  reduced  the  offset 
in  the  LGRAN-C  time  differences.  The  offsets  were  between 
67  and  191  m.  The  drms  was  increased  to  47.3  a  for  the  9940 
X-Y  combination  and  116.9  m  for  the  9940  W-Y  pair. 

The  final  and  most  accurate  method  applies  multiple 
observed  ASF  Correctors  at  five  minute  latitude  and 
longitude  intervals  to  LORAN-C  lines  of  position.  This 
method  again  reduced  the  offset  in  the  time  difference. 
This  offset  was  between  3  and  12  m  for  the  9940  X-Y 
combination.  Part  of  the  offset  may  have  been  a  result  of 
the  microwave  positioning  system.  Reference  is  made  to 
Anderson’s  [  1982  ]  paper  (in  preparation)  on  the  evaluation  of 
multiple  lines  of  position. 


The  drms  values  were  also  reduced  to  38.3  a  for  the  9940 
X-T  rates  and  89.3  a  for  the  9940  H-Y  combination.  The  38.3 
a  1  dras  can  be  decreased  by  improving  the  sampling  time  for 
LORAN-C  receivers.  Nelson  obtained  30  a  1  dras  for  the  9940 
X-Y  rates  with  special  LORAN-C  equipment  used  in  San 
Francisco  Bay,  California  (General  Electric  Co.,  1979]. 
Improving  the  sampling  tiae  for  LOR AN-C  receivers  used  as 
positioning  equipment  for  hydrographic  surveys  should  be 
investigated. 

Hith  dras  values  of  38.3  a  with  the  possibility  of 
obtaining  30.0  a  1  drms  and  offsets  ranging  from  3  to  12  a, 
it  may  be  possible  to  use  LORAN-C  for  hydrographic  surveys 
at  scales  of  1:80,000  or  less  using  multiple  observed  ASF 
Correctors.  The  use  of  DflAHTC  LORAN-C  Correction  Tables 
should  not  be  ignored.  After  updating  these  ASF  Correctors 
with  observed  data,  the  LORAN-C  Correction  Tables  may  allow 
LORAN-C  to  be  used  as  a  positioning  system  for  hydrographic 
surveys. 


i££mil  4 

EIQG  BAfl 


C  PPOGBAH  G*BTLQ 

C 

C  GENE  RAL  PROGRAM  FOR  D  ETER  MINIMS  GP  FROM  KNOWN  STATION 

C  POSITIONS  AID  THE  DISTANCES  FROM  THEM  USING  LEAST 

C  SQUARES.  GP  TO  XT  AND  XT  TO  GP  ARE  DETERMINED  BT 

C  SUBROUTINES  GPTOXT  AID  XTTOGP  WHICH  ARE  BASED  ON  THE 

C  MODIFIED  TRANSVERSE  MERCATOR  PROJECTION  (MTM)  . 

C 

C  PROGRAMMED  BT  GERALD  E.  WHEATON ,  LT.  NOAA 

C 

C  LEAST  SQUARES  ADJUSTMENT  BY  PAUL  R.  WOLF,  PH.D.  AND 

C  REVISED  BT  D.  LEATH,  LCDR 

C 

C  PROGRAM  INPUT  VARIABLE  NAMES 

C  TITLE  *  ANT  JOB  IDENTIFICATION  NAMES  OR  NUMBERS 

C  M  AND  •  THE  NUMBER  OF  EQUATIONS  (M)  AND  UNKNOWNS  (N) 

C  XO  AND  TO  *  BEST  QUESTIMATE  OF  THE  POSTION 

C  STA  (I  r  1)  «  X  COORDINATE  OF  KNOW  STATION 

C  ST A  (1 , 2)  *  T  COORDINATE  OF  KNOW  STATION 

C  STA (1,3)  »  DISTANCE  FROM  KNOW  STATION 

C  A  (I,  J)  »  THE  COEFFICIENT  MATRIX 

C  EL  (I ,  J)  *  THE  CONSTANT  MATRIX 

C  QLL (I , J)  =  THE  WEIGHT  MATRIX  (WEIGHTS  ARE  ENTERED  AS 

C  1'S  IF  THE  SOLUTION  IS  EQUALLY  WEIGHTED) 

C 

IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  XCORD  (30)  ,  Y  CORD  (30)  ,WT  (30)  ,ISNO(30) 

COMMON  /ISTAT/  STA  (30 , 3)  ,  FXY  (30) 

COMMON  /LSQX/  A  (30  ,30  )  ,  EL  ( 30 ,  1)  ,  QLL  ( 30 , 30)  ,  AT  ( 3 0 , 30)  , 
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1AQ  (30,30)  ,  QXX(30f30)  ,  AQL(30,10),  1(30,1),  V(30,1) 
2VAR(30)  ,  TITLE  (80) 

HEAD  AND  WRITE  OUTPUT  TITLE 

SHITE  (6,509) 

READ  (5,710)  TITLE 
SHITE  (6,710)  TITLE 

DEFINE  NUMBER  OF  SIGNALS  (NOT  SHEATEH  THAN  30)  AND 
NUMBER  OF  DATA  SETS. 

READ  (5,  502)  NSIG,NDAT  A 
FORMAT  (15, 1 10) 

DEFINE  THE  CENTRAL  MERIDIAN  (CMER)  ,  FALSE  EASTING 
(FEST),  AND  CENTRAL  LATITUDE  (CLAT)  . 

CMER  IS  EXPRESSED  IN  DEGREES,  MINUTES,  AND  SECONDS. 
FEST  IS  THE  X-COORDINATE  THAT  IS  ASSIGNED  TO  THE 
CENTRAL  MERIDIAN  AND  IS  EXPRESSED  IN  METERS. 

CLAT  IS  DEFINED  AS  THE  CONTROLLING  LATITUDE. 

IT  IS  USED  TO  REFERENCE  THE  Y-COORDIN ATES  AND 
IS  EXPRESSED  IN  METERS. 

READ (5,503)  ILO NC, ILM INC, RLSECC , FEST , CLAT 
FORM AT ( IX, 13, IX, 12, IX ,F8. 5, 1X,F7. 1 ,F10. 1) 

CMER  =  ((IABS  (ILONC)  *  60  ♦  ILMINC)  *  60)  +  RLSECC 

DEFINE  VARIABLE  FOR: 

NUMBER  OF  EQUATIONS  (M) 

NUMBER  OF  UNKNOWNS  (N) 

IPAGE  =  NUMBER  OF  LINES  PER  PAGE. 


IPAGE  =  1 


M  =  4 
N  =  2 
C 

C  NULL  HEIGHTS 

C 

DO  3  IM»1,30,1 

DO  3  JM=1,30,1 

3  QLL(IM,JM)  *  0.0 

C 

C  HEAD  STATION  NUMBERS  , POSITION,  AND  WEIGHTS. 

C  CONVERT  POSITIONS  (GP)  TO  SECONDS  AND  THEN  TO  XY. 

C 

DO  12  J=1 , NSIG  ,  1 

READ  (5,  800)  ISNO  (J)  ,IL  AT  ,1  MIN,  RSEC,  JLON,  JMIN  ,  SSEC,  HT  (J) 
800  F0RMAT(1X,I3,I4,I3,F9 .5,I5,I3,F9. 5,F5. 1) 

C 

RMAST  =  ((lABS(ILAT)  *  60  ♦  IMIN)  *  60)  +  RSEC 
RMASTL  =  ((IABS(JLON)  *  60  ♦  JMIN)  *  60)  +  SSEC 
C 

CALL  GPTOXY  (RMAST,  RMA  STL,  XMETER  ,  YMETER,  FEST  ,  CL  AT  ,  CMER) 
XCORD(J)  =  XMETER 
12  YCORD(J)  »  YMETER 

C 

C  READ  DATA  (STATION  NUMBERS  AND  THE  DISTANCES) 

C 

15  DO  40  JCOUNT=1  ,NDATA  ,  1 

READ  (5, 805)  IF,  IS,  IT,  I4,NREC 
805  FORM  AT  (515) 


DO  16  ICOUNT=1 

,NSI3, 

1 

IF  (ISNO  (I COUNT) 

.EQ. 

IF) 

IF=ICOUNT 

IF  (ISNO  (I COUNT) 

.EQ. 

IS) 

IS=ICOUNT 

IF  (ISNO  (ICOUNT) 

.EQ. 

IT) 

IT=ICOUNT 

IF  (ISNO  (ICOUNT) 

.EQ. 

14) 

I4=IC0UNT 

C 
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STA  (1,1)  *  XCORD  (IF) 

STA(1 ,2)  =  YCORD  (IF) 

QLL(1,1)  =  WT  (IF) 

STA  (  2, 1 )  =  XCORD(IS) 

STA(2f  2)  =  YCORD  (IS) 

QLL(2,2)  =  WT(IS) 

STA  (3,1)  =  XCORD  (IT) 

STA(3,2)  =  YCORD(IT) 

QLL(3,3)  =  WT  (IT) 

STA  (4,1)  =  XCORD  (14) 

STA  (4,2)  =  YCOR  D  (14) 

QLL  (4,4)  =  WT  (1 4) 

READ  THE  DISTANCE  RECORD  AND  LORAN  RATE. 

DO  38  KCOUNT* 1 , NREC, 1 

READ  (5,507)  STA  (1,3)  ,  STA(2,3)  ,STA  (3,3)  ,STA  (4,3) 
1  RATE1,RATE2 
FORM  AT (6 FI  0. 1) 

DETERMINE  BEST  GUESS  COORDINATES  XO  AND  TO 
WITH  SUBROUTINE  GUESS. 

CALL  GUESS  (XO,YO) 

I  JUMP  =  0 

COMPUTE  FXY,  A  AND  L  MATRIX 
DO  25  1=1, M 

FXY  (I)  =  DSQRT  (DABS  (  (  XO-STA  (I,  1)  )  **2 


25 

C 

C 

c 

c 

c 

c 


c 

c 

c 

c 


30 

C 

c 

c 

35 

C 

C 

C 


C 

C 


1  (YO-STA  (I,  2)  )  **2)  ) 

A  {1,1)  *  (XO-ST  A  (1 , 1)  )  /  (FXY(I)) 

A  (I,  2)  =  (YO-STA  (1,2)  )  /  (FXY  (I)  ) 

EL  (1,1)  *  STA  (1 ,3)  -  FXY (X) 

CALL  SUBROUTINE  LSQR 

CALL  LSQR(M,N) 

COMPUTE  THE  NEW  QUESSTIMATE  FOR  XO  AND  YO 

XO  =  XO  ♦  X  (1,1) 

YO  =  YO  +  X  (2 , 1 ) 

EXIT  IF  STANDARD  ARE  MET  USING  IJUMP  OR 
XO  AND  YO  CUT  OFF 

IF  (DABS  (X  (1,  1)  )  .LE.1. 00  .AND.  DABS  (X  (2, 1)  )  .  LE.  1 . 00) 
1  GO  TO  35 

IF  (IJUMP  .EQ.  10)  GO  TO  3  5 
GO  TO  20 

COMPUTE  ERROR  ELLIPSE 

CONTINUE 

CALL  ELIPSE  (SU,  SV) 

CONVERT  XY  TO  GP 

CALL  XYTOGP (XO, YO, SEC LAT, SECLON , FSST , CLAT , CM ER) 

CALL  TODMS  (SECL  AT,  IDE  GP ,1 MINP,  RS  ECP) 

CALL  TODMS  (SECLON,  JDE  GP  ,JMINP ,  S SECP) 

PAGE  AND  CONTINUE  WITH  NEXT  SET  OF  OBSERVATIONS 
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WRITE  (6,505)  IDEGP,IM  INP, RSECP,  JDEGP,  JMINP,  SSECP , 
1  BATE1, BATE2 

WRITE  (6,505)  ID EGP, IMINP, RSECP,  JDEGP, JMINP, SSECP, 
1  SO,SV 

FORM AT (14, 13, F6. 2, 15,  I3,F6 . 2, 2F9. 2) 

IF  (I PAGE  .EQ.  5  0)  WRITE  (6,509) 

IF  (I PAGE  .EQ.  50)  IPAGE  *  0 
IPAGE  =  IPAGE  +  1 

FORMAT  STATEMENTS 

FORMAT  (80  A1) 

FORM  AT(1H1) 

CONTINUE 

STOP 

END 

SUBROUTINE  LSQR  (M,  N) 


IMPLICIT  REAL*8  (A-H,0-Z) 

COMMON  /LSQX/  A  (30  ,30 )  ,  EL  ( 30, 1)  ,  QLL  (30 , 30)  ,  AT  ( 30 , 30) 
1AQ  (3  0,30)  ,  QXX(30,30)  ,  AQL(30,10),  X(30,1),  V(30,1), 
2VAR  ( 30)  ,  TITLE  ( 80) 

COMPUTE  A  TRANSPOSE  BY  TRANSPOSING  THE  A  MATRIX  (AT) 

DO  61  1*1, H 
DO  61  J=1,N 
61  AT  (J,I)  =A  (I,J) 


USING  STEPS  (1),  (2),  AND  (3)  COMPUTE  THE  INVERSE 
OF  THE  TRANSPOSE  (AT)  *  WEIGHTED  MATRIC  (QLL)  * 
MATRIX  A  =  QXX. 


c 

C  (1)  COMPOTE  AQ  =  AT  *  QLL 

C 

DO  71  1*1, N 
DO  7 1  J*1 , M 
AQ  (I,J)  *0. 

DO  71  K*1,M 

71  AQ  (I,J)  *AQ  <I,J)  +  (AT(I  ,K)*QLL(K,  J)  ) 

C 

C  (2)  COMPOTE  QXX  *  AQ  *  A 

C 

DO  81  1*1,  N 
DO  81  J*1,N 
QXX(I,J)=0. 

DO  81  K=1 , M 

81  QXX  (I ,  J)  *QXX  ( I,  J)  +  AQ  ( I,  K)  *A(K,J) 

C 

C  (3)  INVEfiT  QXX  MATH IX 

C 

DO  307  K=1,N 
DO  302  J=1,N 
IF  (J-K)  304,302,  304 
304  QXX(K,J)=QXX  (K,  J) /QXX  (K,  K) 

302  CONTINUE 

QXX(K,K)  =  1./QXX  (K,  K) 

DO  307  1=1, N 
IF  (I-K)  305,307,  305 
3  05  DO  3  03  J=1,N 

IF  (J-K)  306,303,  306 

306  QXX  (I ,J)  =QXX  (I,J) -QXX  (I,K)  *QXX  (K,J) 

303  CONTINUE 

QXX  (I  ,K)  =-QXX  (I  ,  K)  *QX  X  (K,  K) 

307  CONTINUE 
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USING  STEPS  (4)  AND  (5),  COMPUTE  THE  UNKNOWNS  X 
BY  MULT  THE  INVERSE  Q  XX  AND  AQL. 

(4)  COMPUTE  AQL  *  AQ  *  EL 

DO  101  1*1, N 
AQL(I,1) *0. 

DO  101  K=  1 ,  M 

101  AQL  (1,1)  *AQL  (I,  1)  ♦  AQ  ( I,  K)  *EL(K,1) 

(5)  COMPUTE  X  *  QXX  *  AQL 
DO  201  I* 1 , N 

X  (I,  1)  =0. 

DO  201  K=1,N 

201  X  (I,  1)  =X(I,1)  +QXX  (I,  K)  *AQL  (K,  1) 

(6)  COMPUTE  THE  RESIDUAL  (V  *  A  *  X  -EL) 

DO  3  01  1=1,  M 
V(I,  1)=0. 

DO  3  01  K=  1 ,  N 

301  V  (I,1)=V(I,1)  +  A(I,K)*X(K,1) 

DO  1  1*1, M 

1  V  (I,  1)=V  (1,1)  -EL  (1,1) 

COMPUTE  THE  STANDARD  DEVIATION  OF  UNIT  WEIGHT  SIGMA 
DM  -  NUMBER  OF  OBSERVATIONS 
DN  -  NUMBER  OF  KNKNOWNS 

SIGM A=Q . 

DM=M 

DN=N 

DO  382  1=1, M 


382  SIGMA=SIGHA*V  (1,1)  **2*QLL(I,I) 

SIGN A=DSQHT  (SI GMA/ (D M-DN)  ) 

C 

C  COMPOTE  THE  STANDARD  DEVIATION  OF  THE  ADJOSTED  UNKNOWNS 
C  QXX  -  ARE  THE  ELEMENTS  OF  THE  COVARIANCE  MATRIX. 

C 

DO  446  1=1, N 

446  VAR (I)  =DSQRT  (QXX (I , I)  *SIGMA**2) 

C 

510  CONTINUE 
RETURN 
END 

SUBROUTINE  ELIPSE(SU,  SV) 

C - 

C 

C  SOLVE  FOR  THE  SEMI  MAJOR  AND  SEMI  MINOR  AXIS  OF 

C  THE  ERROR  ELLIPSE 

C 

IMPLICIT  REAL*8  (A-H,0-Z) 

COMMON  /LSQX/  A  (30  f 3  0  )  ,  EL  ( 30,  1)  ,  QLL  (30 , 30)  ,  AT  ( 30 , 30 )  , 

1 AQ  (30,30)  ,  QXX(30,30)  ,  AQL(30,10),  X(30,1),  V(30,1), 
2VAR  (30)  ,  TITLE  (80) 

C 

SUS  =  .5*  (QXX  (1,1)  ♦  QXX  (2 , 2)  ♦  DSQRT  (DABS  (QXX  ( 1 ,  1)  - 
1  QXX  (2,2)  +4.0*0 XX  (1,2)  *QXX  (2,  1)  )  )  ) 

SVS  =  .  5*  (QXX  (1,1)  +  QXX  (2,2)  -  DSQRT  (DABS  (QXX  (1  ,  1)  - 
1  QXX  (2,2)  +4. 0*Q  XX  (1,2)*  QXX  (2,1)  )  )  ) 

C 

SU  =  DSQRT  (SUS) 

SV  =  DSQRT  (SVS) 

C 

RETURN 


ii; 


SUBROUTINE  GUESS(XO,YO) 


SUBROUTINE  GUESS  DETERMINES  THE  BEST  GUESS  COORDINATES 
TO  BE  USED  IN  SUBROUTINE  LSQR.  USE  RIGHT  SIDE  RULE  POR 
STATION  ORDER. 

IMPLICIT  REAL*8  (A-H,0-Z) 

COMMON  /ISTAT/  STA  (30 , 3)  , FXY  (30) 

DETERMINE  DISTANCE  BETWEEN  STATION  1  AND  STATION  2 

D=DSQRT  (  (STA  (2,  1) -STA  (1,1)  )  **2+  (STA  (2, 2) -STA  (1,2)  )**2) 

DETERMINE  ANGLE  ALPHA  BETWEEN  XO , YO/STA2/STA1 

ALPH  A=DARCOS  (  (STA  (2,3  )  **2-STA  (1 , 3)  **2+D**2)  / 

1  (2.  0*STA  (2 , 3)  *D) ) 

DETERMINE  ANGLE  BROVO  BETWEEN  X-AXIS  AND  ST A2-STA1 

BR070»DARSIN  (  (STA  (1,2) -STA  (2,2)  ) /D) 

DETERMINE  X  AND  T  LENGTH 

X=STA  (2,3)  *DCOS  (ALPHA  +BROVO) 

Y=ST  A  (2,3)  * DS IN  (ALPHA  +BROVO) 

DETERMINE  XO  AND  YO 


XO  =  ST  A  (2, 1)  +X 
YO  =  STA  (2, 2)  +Y 


ETURN 


SUBROUTINE  GPTOXY(SECLAT,SECLON,XCO,YCO,FEST,CLAT,CMER) 


1 


IMPLICIT  REAL*8  (A-H,0-Z) 

DATA  E2,RKO,A  /. 00676 8658D0,. 99998D0, 6378206. 4D0/ 

DATA  RKGEO, (11, W2  /0.  0  4848 1 368D0 , 0.  1 1422D0, 2 1 . 73607D0/ 
DATA  W3,W4/5 104. 57 33 8 DO ,63  67399.6  89D0/ 

DATA  RADSEC  /.0 000048 48  1368 1 1 1DO/ 

C 

C 

RADLAT  =  SECLAT  *  RADSEC 
SINLAT  =  DSIN  (RADLAT) 

SIN2LA  =  SINLAT  *  SINLAT 
COSLAT  *  DCOS  (RADLAT) 

COS2LA  *  COSLAT  *  COSLAT 
P  =  (CHER  -  SECLON)  /  10000. 0D0 
V  =  A  /  DSQRT  (1  . 0D0  -  E2*SIN2LA) 

TANCON  *  1.0D0  -  SIN2  LA/COS2LA 
S  =  H4  *  (RADLAT  -  SI NL AT* COSLAT/ 1 0. 0D0**6  * 

*  (W3  -CQS2LA*  (W2- Wl*COS2LA)  )  ) 

T1  *  S  *  RKO  -  CLAT 

T2  a  SKGEO  *  CO  SLAT  *  RKO  *  7 

T3  *  12  *  RKGEO  /  2. 0  DO 

T4  *  T3  *  SINLAT 

T5  =  T3  *  RKGEO  *  COS 2LA  /  3. 0D0 

16  *  15  *  TANCON 

T7  *  (4.0D0  ♦  TANCON)  *T5  *RKGEO  *  SINLAT  /  4.0D0 

XCO  *  (T2  +  (T6*P**2))*P  ♦  FEST 

YCO  =  (T7*P**4)  ♦  (T4*P**2)  ♦  T1 

RETURN 

END 
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SUBROUTINE  XYTOGP{XCO  ,YCO, SECLAT, SECLON,FEST,CLATf CMER) 

C - - - 

C 

IMPLICIT  BEAL*8  (A-H,0-Z) 

C 

DATA  E2,A,SR  /O .00676 8658D0 ,6378206. 4D0 , 

1  0.000004  8481 36  8D0/ 

DATA  W1,W2,W3  /O. 2468 2D0, 30. 02335D0, 5078. 64977D0/ 

C 

D  =  CLAT  *  YCO 

WO  =  0.  15704998  1D0/1  0.0D0**6  *  D 
SINWO  =  DSIN(WO) 

COSWO  =  DCOS(WO) 

COS2WO  =  COSWO  *  COSWO 

PHI1  =  WO  ♦  SIN  WO*COS  WO/1  0. 0DQ**6  * 

*  (W3+COS2  WO  *  (W2+W  1*COS2WO)  ) 

PHI2  *  PHI1  /  0 .99998  DO 

PHI3  =  PHI2  /  SR 

Q  =  (XCO  -  FEST)  /  10.0D0**6 

V  =  A  /  DSQRT  (1  -0D0  -  E2*DSIN  (PHI2)  **2) 

T  =  DCOS (PHI2)  *  SR 
C  =  7  *  0.99998  DO 
T 1  =  10 . 0D0**6  /  (T*C ) 

T2  =  (T 1*  1  0.  0D0 **6)  /  (2.  0D0*C) 

T3  =  (T2*10.0D0**6)  /  (3.0D0*C) 

T4  =  (T3*10. 0D0**6)  /  (4. 0D0*C) 

DELLON  =  (T1-Q**2*T3*  (2.0D0*DTAN (PHI2) **2+1 . ODO) ) *Q 
SECLAT  =  ( (  (3.0D0*DTA  N(PHI2)  **2  +  5. 0D0)  *T4)  *Q**2-T2) 

1  *Q**2*DSIN  (PHI2)  ♦  PHI 3 
SECLON  =  CMER  -  DELLON 
RETURN 
END 

SENTRY 


4 
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DATA  SET  EXAMPLE 


PROGRAM  GPBYLQ 


SHIP'S 

LATITUDE 

(D-M-S) 

36  36  40.20 

36  36  42.41 

36  36  44.56 

36  36  46.94 

36  36  49.09 


POSITION 

LONGITUDE 

(D-M-S) 

121  52  48.62 

121  52  48.16 

121  52  47.52 

121  52  46.98 

121  52  46.59 


OBSERVED 

9940-X 

(/usee) 

27508.79 

27508.81 

27508.83 

27508.76 

27508.76 


LORAN  RATES 
9940-1 
Qisec) 

42742.71 

42742.93 

42743.15 

42743. 26 

42743.48 
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jpgsaaii  b 


PROG  SAM  LORAN 

C  PROGRAM  LORAN 

C 

C  PROGRAM  COMPUTES  LINE  OF  POSITIONS  FOR  LORAN-C  USING 

C  THE  SECONDARY  FACTOR  (5F)  BASED  ON  SEAHATER  EM  MODEL 

C 

IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  AXIS  (13)  ,RF  (13) 

DIMENSION  XSLA7  1  (1000)  ,XSLA72  (1000) 

DATA  RHOSEC,PI, UNC07, RN/2. 06264 8062471  DOS, 

1  3. 1415 92 6535 89 8D0 ,299.79  2458D0,1 . 0003 3 8D0/ 

DATA  XMEAN1 , XME AN2 ,7A  R1 ,7  AR2/0. 0  0D0 , 0. 00D0 , 

1  0. 00D0,0 0.0 DO/ 

DATA  XHEAN3,XME AN4/0. 00D0 f0.00D0/ 

C 

CCCCCC****  *********  ********  ************* 

DATA  AXIS/6. 3782064006,6. 378388  D06 , 6. 377397 1 55D06 

1  6.37816  D06, 6. 37816006,6. 3782491 45D06, 6. 3781 65D06 , 

2  6.378166D06,6. 378165  D06 , 6 . 378 1 45D06 ,6. 3775634D06, 

3  6. 37 82 45 DO 6, 6.  378  13  5  0D06/ 

DATA  RF/6. 3565838  006,  2.  97  D02, 2. 991528  1285D02, 

1  2.9825D02,2. 98 247 16 7 427D 02, 2. 93 465D02 , 2 .9825D02 , 

2  2.983  DO 2, 2 . 98  3D0 2, 2.9825  D02,  6. 3562569D06 , 

3  2.983  D02, 2.9826D02/ 

C*****ELLIPSOID  OPTION  NUMBER 

C  1.  CLARKE  1866  8.  MERCURY 

C  2.  INTERNATIONAL  (HAYFORD)  9.  MARSHALL  ISLAND 

C  3.  BESSEL  1841  10.  NA7Y  8D 

C  4.  AND  (AUSTRALIAN)  11.  AIRY 

C  5.  1967  REFERENCE  12.  KRASSOWSKI  1940 


C  6.  CLARKE  1380  ROD  13.  WGS  1972 

C  7.  S AO 
C 

cccccc******* ******** ******  ************* 

C  CC1-2  =  ELLIPSOID  SOMBER  (K) 

C  CC3-5  *  SOMBER  OP  POI STS  ALOSG  COAST  (IREC) 

C 

READ  (5,  100)  K,IREC 
100  FORMAT ( IX, 12 , 14 ) 

TWOPI*2.*PI 

A=AXIS(K) 

F=1./RF  (K) 

IF  (F.LT.3.D-3)  F=(A-1  ./F) /A 
C 

C  READ  MASTER  ASD  SLA7E  STATIOSS  POSITIONS 

C  THE  FIRST  RECORD  IS  THE  BOMBER  OF  MASTER  AND  SLAVE 

C  STATIONS  FOR  THE  PARTICOLAR  CHAIN. 

C 

WRITE  (6,202) 

202  FORMAT  < 1 H 1 ) 

WRITE  (6,201) 

201  FORMAT  (1H  ) 

READ  (5,  105)  I  LA  TM,  I  MI  NM,RSECM,ILONM,  ILMINM  ,  RLSECM 
WRITE (6, 105)  IL ATM, I M INM, RSECM, ILONM, ILMINM, RLSECM 
105  FORM  AT (IX, 13, IX, 12, IX ,F5. 2, IX, 1 4 , 1 X, 12, IX, F5. 2) 

RMAST  =  (  (IABS(ILATM)  *  60  ♦  IMINM)  *  60  +  RSECM)  / 

1  RHOSEC 

IF  (ILATH  .LT.  0)  RMAST  =  -RMAST 

RMASTL  *  (  (IABS  (ILONM)  *  60  ♦  ILMINM)  *60  +RLSECM)  / 
1  RHOSEC 

IF  (ILONM  .GT.  0)  RMASTL  =  TWOPI  -  RMASTL 
C 

C  READ  THE  FIRST  SLAVE  STATION  AND  CHANGE  THE  LATITUDE 

C  AND  LONGITODE  INTO  RADIANS. 


c 

READ  (5,101)  I  LA  TS, IHI  NS , RSECS ,  ILO  NS  ,  ILHINS ,  RLS  ECS  , 

1  DELAT1 

WRITE  (6,101)  IL&TS ,IH INS, RSECS , ILONS, ILHINS, RLSECS, 

1  DEL  AY1 

101  F0RHAT(1X,I3, 1X,I2,1X,F5. 2, IX, 14, IX, 12, IX, F5. 2, F9. 2) 

RSLAV1  *  ( (IAfiS  (ILATS )  *  60  ♦  ININS)  *  60  +  BSECS)  / 

1  RHOSEC 

IF  (ILATS  .LT.  0)  RSLAV1  *  -RSLAV1 

RSLAL1  *  ((IABS (ILONS)  *  60  ♦  ILHINS)  *  60  *  RLSECS)  / 
1  RHOSEC 

IF  (ILONS  .GT.  0)  RSLAL1  *  TWOPI  -  RSLAL1 
C 

C  READ  THE  SECOND  SLAVE  STATION  AND  CHANGE  THE  LATITUDE 

C  AND  LONGITUDE  INTO  RADIANS. 

C 

READ  (5, 101)  I LA TS,IHI NS, R SECS, I LONS, ILHINS, RLS ECS, 

1  DELAY2 

WRITE (6,101)  IL ATS, IM INS, RSECS, ILO NS, ILHINS, RLSECS, 

1  DELAY2 

RSLAV2  »  ( (IABS  (ILATS)  *  60  «•  IHINS)  *  60  «•  RSECS)  / 

1  RHOSEC 

IF  (ILATS  .LT.  0)  RSLAV2  *  -RSLAV2 

RSLAL2  =  (  (IABS  (ILONS)  *  60  +  ILHINS)  *  60  ♦  RLSECS)  / 
1  RHOSEC 

IF  (ILONS  .GT.  0)  RSLAL2  =  TWOPI  -  RSLAL2 

WRITE  (6,201) 

C 

C  READ  THE  RECORDS  POSITIONS  AND  OBSERVED  RATES,  THEN 

C  COHPUTE  THE  TOTAL  DISTANCE. 

C 

DO  5  50  1= 1 , IR EC , 1 

READ  (5,  104)  ILAT,IHIN  ,  RSEC , ILON , ILHIN,  RLSEC, 

1  RATE1 ,  RATE2 
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FORBAT(1X,I3,I3,F6.3,  I5,I3,F6. 3 , 3X, 2F1 0.2) 

CONVERT  LAT  AND  LONS  TO  RADIANS. 

RPOST  =  ((IABS(ILAT)  *  60  ♦  IBIN)  *  60  ♦  RSEC)  / 

1  RHOSEC 

IF  (ILAT  -LT.  0)  RPOST  =  -RPOST 

RPOSTL  =  ((IABS(ILON)  *  60  +  ILBIN)  *  60  +  RLSEC)  / 

1  RHOSEC 

IF  (ILON  .GT.  0)  RPOSTL  =  THOPI  -  RPOSTL 

COBPOTE  DISTANCES  AND  AZIHUTHS  FROB  THE  OBSERVED  POINT 

BASTER 

CALL  INVER1 (A rF , RPOST , RPOSTL, SB  AS T, RBASTL, F AZN , BAZH , 

1  DISTH) 

UDISTH  *  (RN  *  DISTM)  /  ONCOV 
CALL  SECFAC  (ODISTBfSFM) 

SLAVE1 

CALL  INVEH1 (A,F , RPOST , RPO STL, RS L A VI r RSLAL1 , F AZ 1 , BAZ 1 , 

1  DIST1) 

ODIST1  =  (RN  *  DIST1)  /  ONCOV 
CALL  SECFAC  (ODIST  1  ,SF  1) 

SLAVE2 

CALL  INVER1  (A, F, RPOST  ,  RPO  STL,RSLAV2/RSLAL2,FAZ2,BAZ2, 

1  DIST2) 

0DIST2  *  (RN  *  DIST2)  /  ONCOV 
CALL  SECFAC  (ODIST2,SF2) 


C  COMPOTE  THE  FATES  AND  COMPARE  TO  THE  OBSERVED  RATES 
C 

TDM1  =  ODIST1  -  ODIST M  +  SF1  -  SFM  +  DELAY  1 
DIFF1  =  TDM1  -  RAT  El 
C 

TDM2  =  ODIST2  -  ODIST M  +  SF2  -  SFM  +  DELAY2 
DIFF2  *  TDM2  -  RATE2 
C 

C  COMPOTE  THE  LANE  WIDTH  IN  METERS  BASED  ON  EQOATION 

C  4.20  IN  ELECTRONIC  SURVEYING  AND  NAVIGATION  - 

C  LAORILA,  PAGE  94. 

C 


BE  1  = 

DABS (FAZ1 

-  FAZM) 

3R2  = 

DABS (FAZ2 

-  FAZM) 

WIDTH1 

=  (DIFF1 

*  UNCOV  * 

0.5) 

/ 

DSIN  (BR1  * 

0.5) 

WIDTH2 

=  (DIFF2 

*  UNCOV  * 

0.5) 

/ 

DSIN  (BR2  * 

0.5) 

C 

C  WRITE  THE  POSITION  OF  VESSEL,  COMPOTED  RATE,  OBSERVED 

C  RATE,  AND  THE  DIFF  BETWEEN  THEM. 

C 

WRITE  (6,200)  ILAT,  IMI  N,  RSEC,ILON,  ILMIN,  RLSEC,  R  ATE1 , 
DIFF  1  ,  RATE2  ,  DIF  F2 

200  FORMAT (IX, 1 3, 13 ,1X,F6 . 3 ,1 5 , 13 , 1 X, F6. 3, FI  2. 2 , F8 . 2, 

1  F12.2,F8.2) 

WRITE  (6,201) 

C 

C  XMEAN1  AND  XMEAN2  ARE  THE  MEAN  DIFFS  BETWEEN  THE 

C  COMPUTED  RATE  AND  THE  OBERSSRVSD.  XSLAV1  AND  XSLAV2 

C  ARE  THE  STORED  DIFFS. 

C 

XMEAN1  =  XMEAN1  ♦  DIF  FI 
XMEAN2  =  XMEAN2  ♦  DIFF2 
XMEAN3  =  XMEAN3  ♦  WIDTH1 
XMEAN4  =  XMEAN4  ♦  WIDTH2 
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XSLAV  1  (I)  *  DIFP1 
XSLAV2  (I)  =  DIFF2 
550  CONTINUE 
C 

C  COMPUTE  THE  MEAN  AND  STANDARD  DEVIATION 
C 

XMEAN1  *  XMEAN1  /  IREC 

XMEAN2  *  XMEAN2  /  IREC 

XMEAN3  =  XMEAN3  /  IREC 

XMEAN4  =  XMEAN4  /  IREC 

C 

DO  600  I  =  1,IREC,1 

V AR1  a  VAR1  +  (  (XSLAV  1(1)  -  XMEAN1)**2) 

600  VAR2  =  VAR2  ♦  ((XSLAV2(I)  -  XMEAN2)  **2) 

C 

VAR1  =  VAR1  /  (IREC  -  1.0) 

VAR2  =  VAR2  /  (IREC  -  1.0) 

C 

SD 1  =  DSQRT  (VAR  1) 

SD2  =  DSQRT  (VAR  2) 

C 

WRITE  (6,201) 

WRITE  (6,210)  XMEAN1,SD1,XMEAN3 
WRITE  (6,201) 

WRITE  (6, 211)  XMEAN2,S  D2,XMEAN4 
WRITE  (6,202) 

210  FORMAT  ( IX ,  *  SLAVE  #1,  SEAN  =  ',?10.3, 

1  '  STANDARD  DEVIATION  =  ', 

2  F 10 . 3,  ’  DISTANCE  IN  METERS  =  •  ,F10.3) 

211  FORM  AT  ( IX ,  1  SLAVE  #2,  MEAN  =  ',F10.3, 

1  '  STANDARD  DEVIATION  =  *, 

2  F 10 . 3,  *  DISTANCE  IN  METERS  =  ',F10.3) 
STOP 

END 
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SUBROUTINE  SECFAC(UTDIST,  SF) 


THIS  ROUTINE  HI LL  COMPUTE  THE  SEA  SECONDARY  FACTOR 
UTDIST  *  TOTAL  DISTANCE 
SF  =  SECONDARY  FACTOR 


IMPLICIT  RE AL*8  (A-H,  O-Z) 


COEFFICIENTS 


AO  =  129-04323 

A1  *  -0.40758 

A2  =  0.00064576813 


BO  *  2.741282 

B1  =  -0.011402 
B2  =  .00032774815 


IF  (UTDIST  .GT.  537)  SO  TO  10 

SF  =  (  BO  /  UTDIST)  ♦  B1  +  (  B2  *  UTDIST) 

GO  TO  20 

SF  =  (  AO  /  UTDIST)  ♦  A1  +  (  A2  *  UTDIST) 

CONTINUE 

RETURN 


SUBROUTINE  INVER1 (Ar  ? INV, GLAT1 , GLON 1 , GLAT2 , GLON2  ,FAZ, 
1  BAZ  f  S) 


PI 


9 

3 

X 

i 


4 


C  ***  SOLUTION  OF  THE  GEODETIC  INVERSE  PROBLEM  AFTER 
C  ♦**  T.  VINCENT!  MODIFIED  RAINS  FORD'S  METHOD  WITH  HELMERT'S 
C  *♦*  ELLIPTICAL  TERMS.  EFFECTIVE  IN  ANY  AZIMUTH  AND  AT 
C  ***  ANY  DISTANCE  SHORT  OP  ANTIPODAL  STANDPOINT/FOREPOINT 
C  ***  MUST  NOT  BE  THE  GEOGRAPHIC  POLE 
C 

C  ***  A  IS  THE  SEMI-MAJOR  AXIS  OF  THE  REFERENCE  ELLIPSOID 
C  ***  FINV  IS  THE  FLATTENING  (NOT  RECIPROCAL)  OF  THE 
C  ***  REFERNECE  ELLIPSOID  LATITUDES  AND  LONGITUDES  IN 
C  ***  RADIANS  POSITIVE  NORTH  AND  EAST  FORWARD  AZIMUTHS  AT 
c  ***  BOTH  POINTS  RETURNED  IN  RADIANS  FROM  NORTH  GEODESIC 
C  ***  DISTANCE  S  RETURNED  IN  UNITS  OP  SEMI-MAJOR  AXIS  A 
C 

C  ***  PROGRAMMED  FOR  CDC-6600  BY  LCDR  L. PFEIFER  NGS 
C  ***  ROCKVILLE  MD  18FEB75.  MODIFIED  FOR  IBM  SYSTEM  360 
c  ***  BY  JOHN  G  GERGEN  NGS  ROCKVILLE  HD  7507. 

C 

IMPLICIT  RE AL*8  (A-H,0-Z) 

DATA  EPS/0. 5D-1 3/, PI/3. 14 1 5926535898D0/ 

TWOPI=2.*PI 
R=  1.  -FINV 

TU1  =  R*DSIN  (GLAT  1)  /DCOS  (GLAT1) 

7  TU2=R*DSIN  (GLAT  2)  /DCOS  (GLAT2) 

CU 1= 1 ./DSQRT (TU  1*TU1+  1.) 

SU1=CU1*TU1 

CU2=1 ./DSQRT (TU 2*TU2+  1.) 

S  =  CU1*CU2 
BAZ=S*TU2 
FAZ=BAZ*TU1 
X=GLON2-GLON1 
100  SX=D SIN  (X) 


CX=DCOS  (X) 

T01=CU2*SX 

TU2=  B  AZ -S  U 1 *CU2  *CX 
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SY=DSQHT(T01*T01+TO2*T02) 

CY*S*CX+FAZ 

Y=DATAN2(SY,CY) 

SA=S*SX/SY 

C2A=-SA*SA+1. 

CZ=FAZ+FAZ 

IF  (C2A.GT.0.)  CZ *-CZ/C 2A+CY 
E=CZ*CZ*2.-1. 

C=  ( (-3.  *C2A+4  .)  *FI NV+  4.)  *C2A*FINV/16 . 

D=X 

X= ((E*CY»C+CZ) *SY*C+Y) *SA 
X=  (1  .-C)*X*FINV+GL0N2  -GL0N1 
IF  (DABS  (D-X)  .GT.EPS)  GOTO  100 
FAZ*DAT AN2  (TtJ1,TU2) 

IF  (FAZ.GE.  TWOPI)  FAZ* FAZ-T WOPI 
IF  (FAZ.LT.0.D0)  FAZ=F AZ+THOPI 
BAZ=DATAN2(CU1*SX,BAZ*CX-SU1*C02)  +PI 
IF  (B AZ. GE. TWOPI)  3 AZ* BAZ-T WOPI 
IF  (BAZ.  LT.  0.D0)  BAZ=B AZ+T WOPI 
X  =  DSQRT  ((1./S/R-1.)  *C2A+1.)  +  1 . 

X=  (X-2. )  /  X 
C=1.-X 

C=  (X  *X/4.  + 1 . )  /C 
D=  (0 . 375*X*X-  1 .  )  *X 
X= E*CY 
S=  1.  -E-E 

S=  ((  (  (SY*SY*4.-3.)  *S*CZ*D/6.-X)  *D/4.+CZ)  *SY*D+Y)  *C*A*R 

RETURN 

END 

SENTRY 
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D&TA  SET  EXAMPLE  -  P BOG RAM  LORAN 


SHIP'S  POSITION 
LATITUDE  LONGITUDE 

(D-N-S)  (D-N-S) 

36  43  45.800  121  55  27.160 

36  44  3.400  121  55  32.340 

36  44  21.  180  121  55  37.  390 

36  44  37.490  121  55  46.950 

36  44  53.260  121  55  57.710 


SERVED  LORAN  RATES  8  ERRORS  (E) 
9940-y  EY  9940-H  EH 
{jisec)  (jisec)  (jisec)  {jisec) 

42789.34  -0.49  16294.04  -1.06 

42791.13  -0.38  16293.46  -1.10 

42793.04  -0.38  16292.73  -0.99 

42795.13  -0.58  16292.03  -1.06 

42796.93  -0.51  16291.43  -1.27 
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APPENDIX  g 
PROG  HAM  LOPLC 

C  PROGRAM  LOPLC 

C 

C  PROGRAM  COMPOTES  LINE  OP  POSITIONS  FOR  LORAN-C  USING 

C  THE  SECONDARY  FACTOR  (SF)  AND  ADDITIONAL  SECONDARY 

C  FACTOR  (ASF)  THE  SF  AND  ASF  ARE  BASED  ON 
C  SEMI-EMPIRICAL  TD  GRID. 

C 

IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  PHI  (100)  ,  ELON  (100)  ,DISTM  (100)  ,DIST1  (100)  , 

1  DIST2  ( 10  0)  ,  AXIS  (13)  ,  R?(13)  ,AZ1  (100)  ,AZ2(100)  ,AZM  (100) 
DATA  RHOS EC, PI, UNCOV/2. 06 264806247 1D05, 

1  3.141592653589 8D0 ,299. 69 3 00/ 

C 

CCCCCC******* ********  ******  ************* 

DATA  AXIS/6. 378 2064D0 6, 6. 378388  D06, 6. 377397 1 55D06 , 

1  6.37816  D06,6. 378 16D 06,6 . 37824 91 45D06, 6. 3781 65D06 , 

2  6.378166006,6. 378165  D06 ,6. 378 145D06,6. 3775634D06, 

3  6.378245006,6. 378  1350D06/ 

DATA  RF/6. 3565838  D06,  2.97  D02,  2.99  1528  1285D02, 

1  2.9825D02, 2. 98 247 167 427D 02 , 2 . 9 3 465D02 , 2 .9825D02 , 

2  2.983  D02,2. 983D02,2  .9825  D02,  6. 3562569D06 , 

3  2.983  D02, 2.98  26D02/ 

C*****ELLIPS0ID  OPTION  NUMBER 

C  1.  CLARKE  1866  8.  MERCURY 

C  2.  INTERNATIONAL  (HAYFORD)  9.  MARSHALL  ISLAND 

C  3.  BESSEL  1841  10.  NAVY  8D 

C  4.  AND  (AUSTRALIAN)  11.  AIRY 

C  5.  1967  REFERENCE  12.  KRASSOWSKI  1940 

C  6.  CLARKE  1880  MOD  13.  WGS  1972 


C  7.  S AO 
C 

CCCCCC* **************  ******  ************* 
c  CC1-2  *  ELLIPSOID  NUMBER  (K) 

C  CC3-5  *  NUMBER  OF  POINTS  ALONG  COAST  (IREC) 

C 

READ  (5,  100)  K,IREC,ISTA1,  ISTA2 
100  FORHAT(1X,I2,I4  ,213) 

TWOPI=2.*PI 
A-AXIS  (K) 

F= 1. /PF  (K) 

IF  (F.LT.3.D-3)  F*(A-1  ./F) /A 

C 

C  READ  MASTER  STATION  POSITIONS. 

C 

READ  (5,  105)  ILATM, IMINM,RSECM, ILONM, ILMINM,RLSECM 
105  FORMAT(1X,I3, IX ,12 ,1X ,F5.  2, IX, 14 , IX, 12, IX, F5. 2) 

RM AST  »  (  (I ABS  ( ILATM)  *  60  +  ISINH)  *  60  ♦  RSECS)  / 

1  RHOSEC 

IF  (ILATM  .LT.  0)  RMAST  =  -RMAST 

RMASTL  =  ( (IABS  (ILONM)  *  60  +  ILMINM)  *60  ♦RLSECM)  / 
1  RHOSEC 

IF  (ILONM  ,GT.  0)  RMASTL  *  TWOPI  -  RMASTL 
C 

C  READ  THE  FIRST  SLAVE  STATION  AND  CHANGE  THE  LATITUDE 

C  AND  LONGITUDE  INTO  RA DIANS . DELAY  IS  THE  CODING  DELAY 

C  AND  BIAS  IS  THE  OFFSET  IN  MICROSECONDS. 

C 

READ  (5,  101)  ILATS,  IMI  NS  ,R  SECS,  I  LONS,  ILMINS ,  RLSECS  , 

1  DELAY1,3IAS1 

101  FORM  AT ( IX, 13, IX, 12, IX, F5. 2, IX, 1 4, IX, 12, 1X,F5.2 , 

1  F9.  2,F7.3) 

RSLA71  =  ( (IABS  (ILATS )  *  60  +  IHINS)  *  60  ♦  RSECS)  / 
1  RHOSEC 
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IF  (ILATS  .LT.  0)  RSLAV1  *  -RSLAV1 

RSLAL1  *  ( (IABS  (ILOHS )  *  60  ♦  ILHIHS)  *  60  ♦  RLSECS)  / 
1  RHOSEC 

IF  (ILON S  .GT.  0)  RSLAL1  =  THOPI  -  RSLAL1 
C 

C  READ  THE  SECOND  SLAVE  STATION  AND  CHANGE  THE  LATITUDE 
C  AND  LONGITUDE  INTO  RADIANS. 

C 

READ  (5,  101)  ILATS,  IHI  NS ,RSECS, ILONS,  ILMINS ,  RLSECS, 

1  DEL AY2 , BIAS2 

RSLAV2  *  ((IABS  (ILATS)  *  60  ♦  IMINS)  *  60  ♦  RSECS)  / 

1  RHOSEC 

IF  (ILATS  .LT.  0)  RSLAV2  *  -RSLAV2 

RSLAL2  =  ( (IABS  (ILONS)  *  60  +  ILMINS)  *  60  +  RLSECS)  / 
1  RHOSEC 

IF  (ILONS  .GT.  0)  RSLAL2  *  THOPI  -  RSLAL2 
C 

C  READ  COAST  POINT  LAT  AND  LONG  AND  CONVERT  TO  RADIANS. 

C  STORE  LAT  IN  PHI  AND  LONG  IN  ELON. 

C 

DO  5  00  I=1,IRZC,1 

READ  (5,102)  INUM,ILAT,IMIN,3SSC, ILON, ILMIN , RLSEC 
’02  FORM AT ( IX, 14, 14, 13, F7. 3, I 5, 13, F 7. 3) 

PHI  ( I)  =  ((lABS(ILAI)  *  60  +  IMIN)  *  60  ♦  RSEC)  / 

1  RHOSEC 

IF  (ILAT  .LT.  0)  PHI  (I)  =  -PHI  (I) 

ELON  (I)  =  ((IABS(ILON)  *  60  +  ILMIN)  *  60  +  RLSEC)  / 

1  RHOSEC 

IF  (ILON  .GT.  0)  ELON  (I)  =  TWOPI  -  ELON  (I) 

C 

C  COMPUTE  DISTANCE  FROM  COAST  POINT  TO  MASTER  AND 

C  SLAVE  STATIONS.  SET  UP  COMPUTATION  AND  CALL  INVER  1 

C 

PI  =  PHI (I) 
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El  «  ELON(I) 

CALL  INVERl (A, F, Pi, El ,R MAST, RM AST L, AZF, AZB, S) 
DISTM(I)  »  S 
AZM  (I)  «  AZP 

CALL  INVERl (A,F,P1,E1 ,RSLAV1, RSLAL1 , AZF, AZB,  S) 

DIST1  (I)  a  s 
AZ1 (I)  «  AZF 

CALL  INVERl (A, F, PI, El ,RSLAV2, RSLAL2 , AZF,  AZB,  S) 

DIST2  (I)  »  s 
AZ2  (I)  =  AZF 

WRITE  (6,210)  I,INUM,ILAT,IMIN,  RSEC, ILCN, ILMIN,  RLSEC, 
1  AZM  (I)  , AZ1  (I)  ,  AZ2  (I) 

FORM AT (IX, 14, 14 ,14,13 ,F7. 3 , 15 , 13 , F7. 3, 3F10 . 3) 
CONTINUE 

READ  THE  NUMBER  OF  RECORDS  AND  THEIR  POSITIONS. 
COMPUTE  THE  TOTAL  DISTANCE  OVER  LAND  AND  SEA. 

TDISTM  =  TOTAL  DISTANCE  FROM  MASTER  TO  DATA  POINT. 
TDISTS  *  TOTAL  DISTANCE  FROM  SLAVE  TO  DATA  POINT. 


COMPUTE  THE  FORWARD  AND  BACK  AZIMUTHS 

FAZM  AND  FAZS  =  FORWARD  AZIMUTH  TO  THE  MASTER  AND 
SLAVE  STATION. 

EAZM  AND  3AZS  =  BACK  AZIMUTH  TO  THE  MASTER  AND 
SLAVE  STATION. 


READ  (5,  103)  JSEC 
FORM  AT  ( IX  ,  1 4) 


DO  5  50  1=  1 ,  JR  EC  ,  1 

READ (5, 104) ILAT ,IHIN,  RSEC , ILON,  ILMIN , RLSEC , RAT E 1 , RATE2 
FORM  AT (IX, 13, 13 ,F6.3, I5,I3,F6.3,F9.3,F9.3) 


CONVERT  LAT  AND  LONG  TO  RADIANS 


RPOST  ■  ( (IABS  (ILAT)  *  60  ♦  IMIN)  *  60  ♦  RSEC)  / 

1  RHOSEC 

IF  (ILAT  .LT.  0)  RPOST  »  -RPOST 

RPOSTL  *  ((IABS(ILON)  *  60  ♦  ILMIN)  *  60  ♦  RLSEC)  / 

1  RHOSEC 

IF  (ILON  .GT.  0)  RPOSTL  =  THOPI  -  RPOSTL 

COMPUTE  DISTANCES  AND  AZINOTHS  FROM  THE  OBSERVED 
POINT  TO  THE  INTERPOLATED  SHORE  POINT  AND  WRITE. 

ISTA  =  1 

CALL  SHORPT (I REC,RPOS T, RPOSTL, RMAST,RMASTL , UTDISH, 

1  USD 1ST, ULDIST,  FAZ  M,  B  AZM) 

CALL  SECFAC (OTD ISM, USDIST, ULDIST, BAZM,ISTA, SFM) 

CALL  SHORPT (I REC,RPOS T, RPOSTL, RSL AVI , RSLAL 1 , UTDIS 1 , 

1  USD  1ST, ULDIST,  FAZ1,  B  AZ1) 

CALL  SECFAC  (UTD IS  1 , US  DIST  ,  ULDIS T ,  B A Z 1 , 1 STA 1  ,  SF 1 ) 

CALL  SHORPT (I RE C, EPOS T, RPOSTL, RSL AV2 , RSLAL2 , UTDIS2 , 

1  USDIST , ULDIST,  FAZ2,  B  AZ2) 

CALL  SECFAC (UTD IS2 , US DIST , ULDIST , BAZ2 ,ISTA2 , SF2) 

COMPUTE  THE  RATES  AND  COMPARE  TO  THE  OBSERVED  RATES 

TDM1  =  UTDIS1  -  UTDIS M  ♦  SF 1  -  SF M  ♦  DELAY 1  ♦  BIAS1 
DIFF1  =  TDM  1  -  RATE1 

TDM2  *  UTDIS2  -  f’TDIS  M  +  SF2  -  SFM  ♦  DELAY 2  ♦  BIAS2 
DIFF2  =  TDM2  -  RATE2 

COMPUTE  THE  LANE  WIDTH  IN  METERS  BASED  ON  EQUATION 
4.20  IN  ELECTRONIC  SURVEYING  AND  NAVIGATION  -  LAURILA 
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BR1  »  DABS  (FAZ1  -  FAZ  H) 

I F  (B B 1  .GT.  PI)  BB1  =  THOPI  -  BR1 
BR2  »  DABS  (FAZ2  -  FAZ  H) 

IF  (BB2  .GT.  PI)  BB2  =  TWOPI  -  BB2 

WIDTH  1  »  (DIFF1  *  UNCOV  *  0.5)  /  DSIN(BR1  *  0.5) 

WIDTH2  »  (DIFF2  *  UNCOV  *  0.5)  /  DSIN  (BB2  *  0.5) 

WRITE  THE  POSITION  OF  VESSEL,  OBSERVED  RATES, 

AND  THE  DIFFERENCES  BETWEEN  THEM. 

WRITE  (6,200)  ILAT,  IMI  N,  BSEC,  ILON,  ILMIN,  P.LSEC,RATS1 , 

1  DIFF1, RATE2, DI FF2 

FORMAT;  IX, 12, 13, F7. 3, 14 ,1 3,F7. 3 ,  F 10. 2,  F6. 2,  F 10 . 2  ,  F6 . 2) 

XMEAN1  AND  XHEAN2  ARE  THE  MEAN  DIFFS  BETWEEN  THE 
COMPUTED  RATE  AND  THE  OBERSERVED.  XSLAV1  AND  XSLAV2 
ARE  THE  STORED  DIFFS. 

XMEAN1  =  XMEAN1  +  DIF  FI 

XMEAN2  =  XMT AN2  ♦  DIFF2 

XMEAN3  =  XMEAN3  +  WIDTH1 

XMEAN4  a  XMEAN4  ♦  WIDTH2 

XSLAVI(I)  a  DIFF1 
XSLAV2  (I)  =  DIFF2 
CONTINUE 

COMPUTE  THE  MEAN  AND  STANDARD  DEVIATION 

XMEAN1  *  XMEAN1  /  JREC 

X MEAN 2  *  XMEAN2  /  JREC 

XMEAN3  =  XMEAN3  /  JREC 

XMEAN4  =  XMEAN4  /  JREC 


DO  600  I  *  1,JREC,1 

7AR1  *  7AB1  ♦  (  (XSLA7 1  (I)  -  XMEAN1)**2) 

7AR2  =  7AR2  ♦  ((XSLA7  2(I)  -  XMEAN2)  **2) 

7AR1  =  7AR1  /  (JREC  -  1.0) 

7AR2  =  7AR2  /  (JREC  -  1.0) 

SD1  *  DSQRT  (7AR  1) 

SD2  =  DSQRT  (7AR2) 

WRITE  (6  #201) 

WRITE  (6,210)  XM  EAR  1,SD1,XMEAN3 
WRITE  (6,201) 

WRITE  (6,211)  XMEAN2,S  D2,XMEAN4 
WRITE  (6,202) 

FORMAT  (  IX ,  '  SLA7E  #1,  MEAN  =  *,710.3, 

1  •  STANDARD  DE7IATI0N  =  *, 

2  F 10 . 3, *  1  DISTANCE  T.iJ  METERS  *  *,710.3) 
FORM  AT  ( IX  ,  *  SLA7E  *2,  MEAN  =  *,F10.3, 

1  •  STANDARD  DE7IATI0N  =  ', 

2  F 10 . 3,  '  DISTANCE  IN  METERS  =  •  ,?10.3) 
WRITE  (6,201) 

FORMAT  (1H  ) 

FORMAT ( 1 H 1 ) 


SUBROUTINE  SHORPT  (IRE  C,  RPOST  ,  RPOSTL,  RCOllT,  RCONTL,UDIST, 

1  USD,  OLD,  FA,  BA) 
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KNOWN  LATITUDES  AND  LONGITUDES.  THE  SHORE  POINT  IS 
INTERPOLATED  BETWEEN  TWO  KNOWN  POINTS  USING  THE  TOTAL 
DISTANCE  BETWEEN  THE  POSITION  AND  THE  CONTROL  STATION 
AND  THE  AZIMUTH  BETWEEN  THE  SHORE  POINTS  AND  THE 
RECEIVERS  POSITION. 

IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  FAZM(IOO) 

COMMON/SHORE/PHI  (100)  ,  ELO  N  (100)  ,  UNCOV ,RN,  A  ,  F 
DATA  PI/3. 141 59 265358 98D0/ 

TWOPI  =  2.0  *  PI 

CALL  INVER1  (A,F  ,RPOST ,RPOSTL,RCONT,RCONTL,  FA, BA, RDISTT) 
UDIST  =  (RN  *  RDISTT)  /  UNCOV 

RCOMP  =  99999.99 

DO  10  J=1 ,IREC,  1 
PI  =  PHI  (J) 

El  =  ELON  (J) 

CALL  INVER1 (AfF,P1 ,E1 , RCONT, 2C0 NTL , AZF, AZ3, RDISTL) 

CALL  INVER1  ( A ,F , RPOST , RPOSTL, ?1  , El , FAZ, BAZ , RDI3TS) 

USD  =  (RN  *  RDISTS)  /  UNCOV 
ULD  =  (RN  *  RDISTL)  /  UNCOV 
FAZM(J)  =  AZF 

COMPUTE  THE  DIFFERENCE  3ETWEEN  THE  TOTAL  DISTANCE 
(UDIST)  AND  THE  SUMMATION  OF  THE  DISTANCE  OVER  THE  WAT! 
(USD)  AND  THE  DISTANCE  OVER  THE  LAND  (ULD).  IT  THE 
DISTANCE  IS  LESS  THAN  RCOMP,  UPDAT  RCOMP  AND  JSTA. 

JSTA  IS  THE  CLOSEST  POINT  A LON 3  THE  SHORELINE  WHICH  IS 
NEAR  THE  EM  PROPAGATION  PATH. 


FDIFP  =  DABS (ODIST  -  (USD  ♦  OLD)) 

IF  (FDIFF  .GT.  RCOMP)  GO  TO  10 
BCOMP  =  FDIFF 
JSTA  =  J 
CONTINUE 

NOW  DETERHINE  THE  INTERPOLATED  LATITUDE  AND  LONGITUDE 
SHORE  POINT  USING  AZIMUTH  PERCENTAGE. 

IUPPER  =  JSTA  -  1 
ILOWER  =  JSTA  +  1 
AZ MU  =  FAZM (IUPPER) 

AZML  =  FAZM  (ILO HER) 

IF  (AZMU  .GT.  FA  .AND.  FA  .GE.  FAZM  (JSTA))  ICH  =  IUPPER 
IF  (FAZM  (JSTA)  .  GE.  FA  .AND.  FA  .GT.  AZML)  ICH  *  ILOWER 

RADJ  =  1.00  -  DABS  ((FAZM(ICH)  -  FA)  /  (FAZM  (ICH)  - 
1  FAZM  (JSTA))) 

RNWPHI  =  PHI  (JSTA)  ♦  ((PHI  (ICH)  -  PHI  (JSTA)  )  *  RADJ) 

IF  (ELON  (JSTA)  .  LE.  ELON(ICH))  RADJ  =  1.0  -  RADJ 
RNWELN  =  ELON  (JSTA)  ♦  ((ELON  (ICH)  -  ELON  (JSTA))  *  RADJ) 

CALL  TODMS  (RNWPHI,  IDG  ,MIN,  SEC) 

RHOLD  =  TWOPI  -  RNWELN 

CALL  TODMS (RHOLD, IDGL f MIN L, SECL) 

CALL  INVER  1 (A ,F, RNWPHI, RNWELN, RCONT, RCONTL , AZ, BZ , RDISTL ) 
CALL  INVER1 ( A , F ,RPOST , RPO STL, RN WP HI , RNWELN , A Z, BZ , RDISTS ) 
USD  =  (RN  *  RDISTS)  /  UNCOV 
ULD  =  (RN  *  RDISTL)  /  UNCOV 
RETURN 
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THIS  ROOTINE  HILL  COHPOTE  THE  LAND/SEA  SECONDARY  FACTOR 
OTDIST  =  TOTAL  DISTANCE 

OSDIST  =  DISTANCE  OVER  THE  SEA  HATER  PATH 

OLDIST  =  DISTANCE  OVER  THE  LAND  PATH 

AZI  =  AZINOTH  FROM  NORTH. 

RASTER  =  1  (ISTA) 

W  =2 

X  =3 

Y  =4 


IMPLICIT  REAL*8  (A-H,0-Z) 

SI  =  (.795  /  OSDIST)  +  0.439  +  (.  00245  *  OSDIST) 


IF  (OTDIST  .GT.  540)  GO  TO  10 

S2  =  (3.188  /  OTDIST)  -  0.594  +  (.000329  *  OTDIST) 
GO  TO  20 

S2  =  (128.8  /  OTDIST)  +  0.  187  *  (.000652  *  OTDIST) 


IF  (U SDIST  .GT.  540)  GO  TO  30 

S3  =  (3.188  /  OSDIST)  -  0.594  *■  (.000329  *  OSDIST) 
GO  TO  40 

S3  =  (128.8  /  USDIST)  +  0.  187  +  (.000652  *  OSDIST) 


IF  (OLDIST  .GT.  540)  30  TO  50 

S4  =  (3.188  /  OLDIST)  -  0.594  *  (.000329  *  OLDIST) 
GO  TO  60 

S4  =  (128.8  /  OLDIST)  +  0.187  ♦  (.000652  *  OLDIST) 


5 5  =  1.  428  ♦  (.  00158  *  OTDIST) 

56  =  1.  428  ♦  (.  00158  *  OLDIST) 


TAZI  *  2-0  *  AZI 

GOTO  (70,80,90,100)  ,  ISTA 

70  RHOLD*  (1.010*DSIN  (AZI)  )  -  ( .  1 96*DCOS  (AZI) ) 

1  -  (. 893 *DSIN (TAZI)  )  -  ( . 355*DCOS  (TAZI)  ) 

GO  TO  200 

80  BHOLD  ~  (.323*DCOS  (AZI)  )  -  (.  71 1  *DSIN  (TAZI) ) 

GO  TO  200 

90  RHOLD  =  (.942*DCOS  (TAZI)  ) 

GO  TO  200 

100  RHOLD  =  (.588*DSIN  (TAZI)  ) 

C 

200  S5  =  S5  +  RHOLD 

S6  *  S6  ♦  RHOLD 

SF  =  0.5  *  (S5  ♦  S6  -  SI  +  S2  ♦  S3  -  S4) 

WRITE  (6,500)  SI ,S2,S3  ,S4, S4,S5, S6 ,SF 
500  FORHAT(1X,7F15.  5) 

RETURN 

END 

SUBROUTINE  INVER1  (A,  FINV,  GLAT1 ,  GLON 1  ,GLAT2 , GLON2  ,FAZ, 

1  BAZ  ,  S) 

C - 

See  Appendix  B  for  subroutine  INVERT. 
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r.-: 


b 


Short  points  around  Monterey  Bay,  California.  The  poin 
are  used  to  interpolate  geodetic  point  for  computation 
of  the  distance  oyer  land  and  the  distance  over  sea. 


NO.  LATITUDE  LONGITUDE 


1 

36 

57 

18, 

.606 

122 

05 

37. 

525 

2 

36 

56 

59. 

.264 

122 

03 

01. 

817 

3 

36 

57 

05. 

.076 

122 

01 

31. 

701 

4 

36 

57 

49. 

,538 

122 

01 

07. 

857 

5 

36 

57 

17. 

.949 

121 

58 

19. 

830 

6 

36 

58 

08. 

.589 

121 

57 

07. 

288 

7 

36 

58 

32. 

.140 

121 

55 

10. 

083 

8 

36 

58 

01. 

.498 

121 

53 

57. 

390 

9 

36 

56 

46. 

.115 

121 

52 

22. 

313 

10 

36 

55 

38. 

.140 

121 

51 

24. 

399 

11 

36 

53 

13. 

.806 

121 

49 

46. 

743 

12 

36 

49 

38. 

.384 

121 

47 

48. 

8  95 

13 

36 

47 

39. 

.241 

121 

47 

10. 

818 

14 

36 

46 

27. 

.554 

121 

47 

39. 

6  37 

15 

36 

44 

56. 

.717 

121 

47 

52. 

,416 

16 

36 

41 

14. 

.439 

121 

48 

32. 

642 

17 

36 

39 

17. 

.211 

121 

49 

28. 

533 

18 

36 

37 

31. 

.128 

121 

50 

31. 

728 

19 

36 

36 

23. 

.446 

121 

51 

34. 

833 

20 

36 

36 

03. 

.628 

121 

52 

50. 

879 

21 

36 

36 

24, 

.782 

121 

53 

48. 

453 

22 

36 

37 

18. 

.151 

121 

54 

11. 

628 

23 

36 

38 

00, 

.300 

121 

55 

57. 

538 

4 


4 


4 
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DATA  SET  EXAMPLE  -  PROGRAM  LOPLC 


Jj 

i 


i 

►  _ 

», . 

s 

k  _■ 

» 

i 

w'  • 

4 


SHIP'S 

POSITION 

OBSERVED  LORAN  RATES  &  ERRORS  (E) 

LATITUDE 

(D-M-S) 

LONGITUDE 

(D-M-S) 

9940-Y 

(jisec) 

ET 

(/isec) 

9940-W 
(p  sec) 

EH 

Qusec) 

36 

43 

45.800 

121 

55 

27. 160 

42789.34 

1.  14 

16294.04 

1.  32 

36 

44 

3.400 

121 

55 

32.340 

42791.13 

1.25 

16293.46 

1.29 

36 

44 

21.  180 

121 

55 

37.390 

42793.04 

1.25 

16292.73 

1.40 

36 

44 

37.490 

121 

55 

46.950 

42795. 13 

1.05 

16292.03 

1.33 

36 

44 

53.260 

121 

55 

57.710 

42796.93 

1 . 12 

16291.43 

1.  11 

p 


Y 


« 
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APPENDIX  D 


PROGRAM  LORTAB 


C  PROGRAM  LORTAB 

C 

C  PROGRAM  COMPOTES  LINE  OP  POSITION  FOR  LORAN-C  USING 

C  SF  SALT  WATER  CORRECTION  PACTOR  AND  DHAHTC  CALCULATED 

C  OR  FIELD  OBSERVED  ASP  CORRECTIONS. 

C 

C 

IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  AXIS  (13)  ,RF  (13) 

DIMENSION  XSLAV1  (1000)  ,XSLAV2  (1000) 

DATA  RHOSEC,PI, UNCOV,  RN/2 . 06264 8062471 D05 , 

1  3. 1415 92 6535  898D0 ,29  9.  792458D0 , 1 . 000338D0/ 

DATA  XMEAN1,XMEAN2,VAR1,VAR2/0. 00D0,0.00D0, 

1  0.00D0r0 0.0D0/ 

DATA  XMEAN3,XMEAN4/0. 00D0,0.00D0/ 

C 

CCCCCC**  ***********  ********  ************* 

DATA  AXIS/6. 3782064D06, 6. 378388  D06 , 6. 3773971 55D06 , 

1  6.37816  D06,6. 378 16D 06,6 . 3782U9145D06,6. 378165D06, 

2  6.378166006,6.  378165  D06 , 6. 378 1 4 5D06r 6 . 3775634D06, 

3  6.378245006,6. 378135  0DO6/ 

DATA  RF/6. 3565838  D06,  2.97  D02, 2 .99 1528 12  85D02  , 

1  2.98  25D0  2, 2.  98 247 167 427D 02, 2. 93465D02 ,2 .9825D0 2, 

2  2.983  D02, 2. 983002,2.9825  D02, 6. 3562569D06 , 

3  2.983  D02, 2. 98  26D02/ 

C*****ELLI?SOID  OPTION  NUMBER 

C  1.  CLARKE  1866  8.  MERCURY 

C  2.  INTERNATIONAL  (HAYP7RD)  9.  MARSHALL  ISLAND 

C  3.  BESSEL  1841  10.  NAVY  8D 


c 

4. 

AND  (AUSTRALIAN) 

11.  AIRY 

c 

5. 

1967  REFERENCE 

12.  KRASSOWSKI  1940 

c 

6. 

CLARKE  1880  MOD 

13.  WGS  1972 

c 

7. 

S  AO 

c 

CCCCCC*******  ******  ********  **414141****4141*41 


C  CC1-2  *  ELLIPSOID  NUMBER  (K) 

C  CC3-5  *  NUMBER  OF  POINTS  ALONG  COAST  (IREC) 

C 

READ  (5,  100)  K,IREC,ID1,ID2 
100  FORMAT  (IX, 12, 14, 213) 

TWOPI=2.*PI 
A=  AXIS  (K) 

F=1./RF  (K) 

IF  (F.LT.3.D-3)  F=(A-1./F)/A 
C 

C  READ  MASTER  AND  SLAVE  STATIONS  POSITIONS 

C  THE  FIRST  RECORD  IS  THE  NUMBER  OF  MASTER  AND  SLAVE 

C  STATIONS  FOR  THE  PARTICULAR  CHAIN. 

C 


WRITS  (6,202) 

202  FORMAT  (1H1) 

WRITE  (6,201). 

201  FORMAT  (1H  ) 

READ  (5,  105)  ILATM,  IMI  NM,RSECM,  ILONM,  ILMINM ,  RLSECM 
WRITE (6,105)  IL ATM, IMINM,RSECM, ILONM, ILMINM, RLSECM 
105  FORM  AT (IX, 13, IX, 12, IX, F5.  2, IX, 1 4, IX, 12, IX,  F5. 2) 

RMAST  =  (  (IABS  (ILATM)  *  60  ♦  IMINM)  *  60  +  RSECM)  / 

1  RHOSEC 

IF  (ILATM  .LT.  0)  RMAST  =  -RMAST 

RMASTL  =  ((IABS  (ILONM)  *  60  +  ILMINM)  *60  ♦ RLSECM)  / 
1  RHOSEC 

IF  (ILONM  .GT.  0)  RMASTL  =  TWOPI  -  RMASTL 


BEAD  THE  FIRST  SLAVE  STATION  AND  CHANGE  THE  LATITUDE 
AND  LONGITUDE  INTO  RADIANS. 


READ  (5,  101)  ILATS, INI  NS, RSECS, ILONS, ILMINS,  RLSECS, 

1  DELAY1 

WRITE  (6,1 01)  ILATS ,IH INS, R SECS, ILONS , ILMINS, 3LS ECS, 

1  DEL AY1 

F0BHAT(1X,I3,1X,I2,1X,F5. 2, IX, 1 4, IX, 12, IX, F5. 2, F9. 2) 
RSLAV1  =  ( (IABS  (ILATS)  ♦  60  ♦  IMINS)  *  60  ♦  RSECS)  / 

1  RHOSEC 

IF  (ILATS  .LT.  0)  RSLAV1  *  -RSLAV1 

RSLAL1  =  ((IABS  (ILONS)  *  60  ♦  ILMINS)  *  60  ♦  RLSECS)  / 
1  RHOSEC 

IF  (ILONS  .GT.  0)  RSLAL1  »  TWOPI  -  RSLAL1 


READ  THE  SECOND  SLAVE  STATION  AND  CHANGE  THE  LATITUDE 
AND  LONGITUDE  INTO  RADIANS. 


READ  (5,  101)  ILATS,  IMINS, RSECS,  I  LONS,  ILMINS,  RLSECS, 

1  DELAY2 

WRITE  (6,101)  ILATS ,IM INS, RSECS, ILONS , ILMINS , RLSECS, 

1  DELAY2 

RSLA  V2  =  ( (IABS  (ILATS)  *  60  ♦  IMINS)  *  60  +  RSECS)  / 

1  RHOSEC 

IF  (ILATS  .LT.  0)  RSLAV2  *  -RSLAV2 

RSLAL2  *  ( (IABS  (ILONS)  *  60  +  ILMINS)  *  60  +  RLSECS)  / 
1  RHOSEC 

IF  (ILONS  .GT.  0)  RSLAL2  =  TWOPI  -  RSLAL2 
WRITE  (6,201) 


READ  THE  RECORDS  POSITIONS  AND  OBSERVED  RATES,  THEN 
COMPUTE  THE  TOTAL  DISTANCE. 


DO  550  1=  1 , 1 3 EC  ,  1 


K 


READ  (5, 104)  TLAT,IMIN , RSEC, ILON , ILMIN , RLSEC , 

1  RATE1 ,  RATE2 

104  F0RMAT(1X,I3,I3,F6.3, 15 ,1 3 ,F6 . 3 , 3X, 2F10. 2) 

C 

C  CONVERT  LAT  AND  LONG  TO  RADIANS. 

C 

RPOST  =  ((IABS(ILAT)  *  60  ♦  IMIN)  *  60  ♦  RSEC)  / 

1  RHOSEC 

IF  (ILAT  .LT.  0)  RPOST  =  -RPOST 

RPOSTL  *  ( (I ABS  (ILON)  *  60  *  ILMIN)  *  60  +  RLSEC)  / 

1  RHOSEC 

IF  (ILON  .GT.  0)  RPOSTL  =  TWOPI  -  RPOSTL 
C 

C  COMPOTE  DISTANCES  AND  AZIMUTHS  FROM  THE  OBSERVED  POINT 

C 

C  MASTER 

C 

CALL  INVER1 (A,F , RPOST , RPOSTL, RMAST,RMASTL, FAZM,3AZM, 

1  DISTM) 

ODISTM  =  (RN  *  DISTM)  /  ONCOV 
CALL  SSCF AC  (ODISTM, SF  M) 

C 

C  SLAVE1 

C 

CALL  INVER1 (A,f , RPOST , RPO STL , RS LA V 1 , RSLAL 1 , FAZ 1 , BAZ 1 , 

1  DIST1) 

0DIST1  =  (RN  *  DIST1)  /  ONCOV 
CALL  SECFAC  (ODIST1  ,SF  1) 

C 

C  SLAVE2 

C 

CALL  INVER1 ( A, F , RPOST , RPO STL, RS L A V2 , RSL AL2 , F AZ2 , 3AZ2 , 

1  DIST2) 

0DIST2  =  (RN  *  DIST2)  /  UNC07 
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CALL  SECFAC  (0DIST2,SF  2) 


C 

C  DETERMI NE  THE  ADDITIONAL  SECONDARY  CORRECTORS  FROM  THE 

C  LORAN-C  CORRECTION  TABLE  FOR  THE  WEST  COAST  CHAIN  9940 

C 

CALL  TABLE  (RPOST,RPOS  TL,ID1,ASF1) 

CALL  TABLE (R?OST, RPOS  TL ,ID2, ASF2) 

C 

C  COMPOTE  THE  RATES  AND  COMPARE  TO  THE  OBSERVED  RATES 

C 

TDM1  =  ODIST1  -  ODIST M  +  SF1  -  SFM  ♦  ASF1  +  DELAY1 
DIFF1  =  TDM  1  -  RAT El 
C 

TDM2  =  ODIST2  -  ODISTM  +  SF2  -  SFM  +  ASF2  ♦  DELAY2 
DIFF2  =  TDM2  -  RATE2 
C 
C 

C  COMPOTE  THE  LANE  WIDTH  IN  METERS  BASED  ON  EQOATION 

C  4.20  IN  ELECTRONIC  SORVEYING  AND  NAVIGATION  - 

C  LAORILA,  PAGE  94. 

C 

BE  1  =  DABS  (FAZ1  -  FAZM) 

BR2  =  DABS  (FAZ2  -  FAZM) 

WIDTH1  =  (DIFF1  *  ONCOV  *  0.5)  /  DSIN(BR1  *  0.5) 

WIDTH2  =  (DIFF2  *  UNCOV  *  0.5)  /  DSIN(BR2  *  0.5) 

C 

C  WRITE  THE  POSITION  OF  VESSEL,  COMPOTED  RATE,  OBSERVED 

C  RATE,  AND  THE  DIFF  BETWEEN  THEM. 

C 

WRITE  (6,200)  IL  AT,  IMI  N,  RSEC,ILON,  ILMIN,  RLSEC,RAT31 , 
DIFF 1 , RATE2 ,DIF  F2 

200  FORM  AT (IX, 13, 13 ,1X,F6 . 3 ,1 5 ,13, 1 X, F6. 3 ,F12. 2 , F8 . 2 , 

1  F12.2,F8.2) 

WRITE  (6,201) 
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XMEAN1  AND  XMEAN2  ARE  THE  MEAN  DIFPS  BETWEEN  THE 
COMPUTED  RATE  AND  THE  OBERSERVED.  XSLAV1  AND  XSLAV2 
ARE  THE  STORED  DIFFS. 


XMEAN1  =  XMEAN1  +  DIF  FI 
XMEAN2  a  XMEAN2  ♦  DIFF2 
XMEAN3  =  XMEAN3  ♦  BIDTH1 
XMEAN4  =  XMEAN4  ♦  WIDTH2 
XSLAV  1  (I)  =  DIFF1 
XSLAT2  (I)  =  DIFF2 
CONTINUE 


COMPUTE  THE  MEAN  AND  STANDARD  DEVIATION 


XMEAN1  = 
XMEAN2  = 
XMEAN3  = 
XMEAN4  = 


XMEAN1  /  IREC 
XMEAN2  /  IREC 
XMEAN3  /  IREC 
XMEAN4  /  IREC 


DO  600  I  =  1,IREC,1 

VAR1  =  VAR1  ♦  ((XSLAV  1(1)  -  XMEAB1)**2) 

VAR2  =  VAR2  ♦  ((XSLAV  2  (I)  -  XMEAN2)**2) 

VAR1  =  VAR1  /  (IREC  -  1.0) 

VAR2  =  VAR2  /  (IREC  -  1.0) 

SD1  »  DSQRT  (V  AR  1) 

SD2  =  DSQRT  (VAR2) 


WRITE  (6,201) 

WRIT  E  (6 , 2 1  0)  XMEAN  1,SD1,XMEAN3 
WRITE  (6,201) 

WRITE  (6 , 2 1  1 )  XMEAN2,SD2,XMEAN4 


:»] 


WRITE  (6,202) 

210  FORMAT  ( IX,  '  SLAVE  #1,  MEAN  =  »,F10.3, 

1  «  STANDARD  DEVIATION  =  ', 

2  F10.3,  »  DISTANCE  IN  METERS  =  *,F10.3) 

211  FORMAT  ( IX, '  SLAVE  #2,  MEAN  =  »,F10.3, 

1  •  STANDARD  DEVIATION  =  •  , 

2  F 10 . 3,  ’  DISTANCE  IN  METERS  =  •  ,F10.3) 
STOP 


C 

C 

c 

c 

c 

c 

c 

c 

c 


c 


c 


10 

c 


SUBROUTINE  SECF  AC ( UTD  1ST,  SF) 


THIS  ROUTINE  WILL  COMPUTE  THE  SEA  SECONDARY  FACTOR 
UTDIST  =  TOTAL  DISTANCE 
SF  =  SECONDARY  FACTOR 

IMPLICIT  RE AL*8  (A-H,0-Z) 


COEFFICIENTS 

AO  =  129.04323 

A 1  =  -0.40758 

A2  =  0.00064576813 

BO  =  2.741282 

B 1  =  -0.011402 
B2  =  .000327748  15 

IF  (UTDIST  .GT.  537)  GO  TO  10 

SF  =  (  BO  /  UTDIST)  +  B1  ♦  (  B2  *  UTDIST) 

GO  TO  20 

SF  =  (  AO  /  UTDIST)  «■  A1  +  (  A2  *  UTDIST) 
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SUBROUTINE  TABLE  (RL&T  ,RLON,ID,  ASF) 


SUBROUTINE  TABLE  SELECTS  THE  PROPER  ASF  CORRECTOR  FROM 
THE  LORAN-C  CORRECTION  TABLE  PUBLISHED  BY  THE  DEFENSE 
MAPPING  AGENCY. 

RLAT  =  POSITION  LATITUDE  IN  SECONDS 
RLON  *  POSITION  LONGITUDE  IN  SECONDS 
ID  *  LORAN-C  CHAIN  IDENTIFIER 
a  =  i 
x  =  2 

Y  =  3 

ASF  a  ADDITIONAL  SECONDARY  FACTORS 


THE  FOLLOWING  TABLES  OF  ASF  CORRECTORS  ARE  FOR 
MONTEREY  BAY,  CALIFORNIA  -  9940  -W,  -X,  -Y. 

IMPLICIT  REAL*8  (A-H,  O-Z) 

DIMENSION  TABLE  W  (3  ,5)  ,TAB  LEX  (3 , 5)  ,  T ABLEY  (3 , 5) 


DATA  TABLEW/  1. 6D0 , 1 . 6D0 , 0. 0D0 , 

1  1.  600,1.  500,1.400, 

2  1.  500,1.400,1.600, 

5  1. 400,1.  3D0,1 .500, 

»  1. 3D0 ,0.  ODO, 0.0D0/ 

DATA  TABLEX/  -0 . 9D0, - 0. 9D0 ,  0.0D0, 
-0.9D0,- 0.9D0,-0.8D0, 
-1  .ODO,-  1. ODO, -1. ODO, 
-1 .ODO,-  1.200,-1.100, 


1 

2 

3 
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1  .  IDO#  O.ODO,  O.OD 0/ 


DATA  TABLET/  0. 2D0 , 0. 2DO, 0 . ODO, 

1  0-  3D0,0.  3DO,0.3D0, 

2  0.  3D0,0.  2D0,0.4D0, 

3  0.  4 DO  ,0.  3D0,0.6D0, 

4  0. 5D0f0. ODO, O.ODO/ 

DATA  RH0SEC,PI/2. 0626 4806 2471035,3.1 4 159 265 358 98 DO/ 

TWOPI  *  PI  *  2.0 

CONVEB  SLAT  AND  RLON  TO  SECONDS 

HLAT  =  SLAT  *  8H0SEC 
HLON  =  TWOPI  -  SLON 
HLON  =  HLON  *  RHOSEC 

STARTING  LAT  AND  LONG  POE  SEASCH 

LAT  »  37/00/00.0  LONG  =  122/05/00.0 


DETERMINE  THE  ASF  COS  RECTOR  FOR  THE  LORAN-C  COMBINATION 


LATITUDE 

SLAT  =  133200.0 
SLON  =  439500.0 
RMID  =  300.0 
RDIFF  =  150.0 

DO  10  J= 1 , 5 , 1 
SLAT  =  SLAT  -  RMID 
ULAT  =  SLAT  ♦  RDIFF 
VLAT  =  SLAT  -  RDIFF 


IF  (hLAT  .LT.  OLAT  .AND.  hLAT 
CONTIN0S 


-3  9 


FLAT)  GO  TO  15 


LONG  I  TO  DE 

DO  30  1=1,3, 1 
SLON  =  SLON  -  RHID 
OLON  =  SLON  ♦  RDIFF 
FLON  =  SLON  -  RDIFF 

30  IF  (hLON  .LT.  OLON  .AND.  h  LON  .3T.  FLON)  GO  TO  35 
35  CONTINUE 

DETERMINE  ASF  CORRECTOR 

ASF  =0.0 

IF  (ID  .EQ.  1)  ASF  =  T  ABLER  (I,  J) 

IF  (ID  .EQ.  2)  ASF  =  TABLEX(I,J) 

IF  (ID  .EQ.  3)  ASF  =  TABLEY  (I,J) 

RETORN 
END 

SOBROOTINE  INFER1 ( A, F INF, GLAT 1 r GLON 1 , GLAT2 r GLON2 ,FAZ , 


DATA  SET  EXAMPLE  -  PROGRAM  LORTAB 


SHIP'S  POSITION 
LATITUDE  LONGITUDE 

(D-M-S)  (D-M-S) 


36  43  45.800  121  55  27.  160 


36  44  3.400  121  55  32.  340 


OBSERVED  LORAN  RATES  &  ERRORS  (E) 
9940- Y  El  9940-W  EW 

(psec)  (ji  sec)  (jisec)  (jisec) 

42789.34  -0.29  16294.04  0.34 

42791.13  -0.18  16293.46  0.30 

42793.04  -0.  18  16292.  73  0.41 

42795.13  -0.38  16292.03  0.34 

42796.93  -0.31  16291.43  0.  13 


APPENDIX  E 
PROGRAM  ASFSBL 


C  PROGRAM  ASFSEL 

C 

C  PROGRAM  DETERMINES  OBSERVED  ASF  CORRECTORS  BY  SCANNING 

C  DATA  AT  1  DEGREE  LAIITUDE  AND  LONGITUDE  INTERVALS.  THE 

C  ASF  CORRECTORS  ARE  DETERMINED  BY  SUBTRACTING  THE 

C  CALCULATED  TD  USING  THE  SEAWATER  SECONDARY  FACTOR  FROM 

C  THE  OBSERVED  TD  RATES. 

C 

IMPLICIT  REAL *8  (A-H,0-Z) 

DIMENSION  AXIS  ( 13)  #RF  (13) 

DIMENSION  ASFCR1  (16,26)  ,  ASFCR2  ( 16  ,26)  , 

1  INOI  (16,26)  ,IN02(16, 26) 

DATA  RHOSEC, PI, UNCOV, RN/2 . 06264 8062471D05, 

1  3. 1415 92 6535 89  8D0 ,29  9. 79  2458D0 , 1 . 000338D0/ 

DATA  XMEAN 1 , XME AN2 , V A  R1 ,V AR2/0. 00D0,0.00D0, 

1  0.00D0,00.0D0/ 

DATA  XHEAN3 ,XME AN4/0.  00D0,0.00D0/ 

C 

CCCCCC* **************  ******  ************* 

DATA  AXIS/6. 3782064D0  6,6.  378388  D06, 6. 3773971 55D06 , 

1  6.37816  D06, 6. 37816D06, 6. 378249145006,6. 378165D06, 

2  6.378166D06,6. 378165  D06 , 6 . 378 1 45D06 ,6 . 377 5634D06 , 

3  6.378245006,6.  3781350D06/ 

DATA  RF/6. 3565838  D06,2.97  D 02 , 2 . 99 152 8 12 85D02 , 

1  2.98  25D0  2,2.98  247167  427D02,2.93465D02,2.9  8  25D02, 

2  2.983  D02,2. 983D02,2 .9825  D02, 6. 3562569D06 , 

3  2.983  D02,  2.98  26D02/ 

C*****ELLI?SOID  OPTION  NUMBER 

C  1.  CLARKE  1866  8.  MERCURY 


C  2.  INTERNATIONAL  (HAT FORD)  9.  MARSHALL  ISLAND 
C  3.  BESSEL  1841  10.  NAVY  8D 

C  4.  AND  (AUSTRAL IAN)  11.  AIRY 

C  5.  1967  REFERENCE  12.  KRASSOBSKI  1940 

C  6.  CLARKE  1880  BOD  13.  WGS  1972 

C  7.  S AO 

C 

CCCCCC* ************** ******  ************* 

C  CC 1-2  =  ELLIPSOID  NUMBER  (K) 

C  CC3-5  *  NUMBER  OF  POINTS  ALONG  COAST  (IREC) 

C 

READ  (5,  100)  K,IREC,ID1,ID2 
100  FORMAT(1X, 12,14,213) 

TWOPI=2.*PI 
A= AXIS  (K) 

F=1./RF  (K) 

IF  (F.LT.3.D-3)  F*(A-1./F)/A 
C 

C  READ  MASTER  AND  SLAVE  STATIONS  POSITIONS 

C  THE  FIRST  RECORD  IS  THE  NUMBER  OF  MASTER  AND  SLAVE 

C  STATIONS  FOR  THE  PARTICULAR  CHAIN. 

C 

WRITE  (6,202) 

202  FORMAT  ( 1 H 1 ) 

WRITE  (6,201) 

201  FORMAT (1H  ) 

READ  (5,  105)  I  LATH,  IMI  NM,RSECM,ILONM,  ILMINM  ,  RLSECM 
WRITE  (6,105)  IL ATM ,IM INM, RSECM, ILONM, ILMINM , RLS ECM 
105  FORM  AT (IX, 13, IX, 12, IX ,F5. 2, IX, 14, IX, 12, IX, F5.2) 

RMAST  *  (  (IABS  (ILATM)  *  60  ♦  IMINM)  *  60  ♦  RSECM)  / 

1  RHOSEC 

IF  (ILATM  .LT.  0)  RMAST  =  -RMAST 

RMASTL  *  (  (IABS  (ILONM)  *  60  ♦  ILMINM)  *60  ♦RLSECM)  / 
1  RHOSEC 


IF  (ILONM  .GT.  0)  8 MAS TL  =  TWOPI  -  BMASTL 
C 

C  READ  THE  FIRST  SLAVE  STATION  AND  CHANGE  THE  LATITUDE 

C  AMD  LONGITUDE  INTO  RADIANS. 

C 

READ  (5,  101)  ILATS,IMINS,RSECS,ILONS,ILMINS,RLSECS, 

1  DEL ATI 

WRITE  (6,101)  IL  ATS  ,IM INS  , RSECS,  I  LONS  ,  ILMINS  ,  RLSECS  , 

1  DELAT1 

101  FORM AT (IX, 13, 11,12, IX, F5. 2, IX, 1 4 , IX, 12, IX, F5 . 2, F9 . 2) 

RSLAV1  *  ( (IABS  (ILATS )  *  60  ♦  IMINS)  *  60  ♦  RSECS)  / 

1  RHOSEC 

IF  (ILATS  .LT.  0)  RSLAV1  *  -RSLAV1 

RSLAL1  »  ((IABS  (ILONS)  *  60  ♦  ILMINS)  *  60  ♦  RLSECS)  / 
1  RHOSEC 

IF  (ILONS  .GT.  0)  RSLAL1  =  TWOPI  -  RSLAL1 
C 

C  READ  THE  SECOND  SLAVE  STATION  AND  CHANGE  THE  LATITUDE 

C  AND  LONGITUDE  INTO  RADIANS. 

C 

READ  (5,  101)  ILATS, IMINS, RSECS, ILONS, ILMINS, RLSECS, 

1  DEL  AT2 

WRITE  (6,101)  ILATS, I  SINS,  RSECS,  ILONS,  ILMINS,  RLSECS, 

1  DELAT2 

RSLAV2  *  (  (IABS  (ILATS)  *  60  +  IMINS)  *  60  +  RSECS)  / 

1  RHOSEC 

IF  (ILATS  .LT.  0)  RSLAV2  =  -RSLAV2 

RSLAL2  =  ( (IABS  (ILONS)  *  60  ♦  ILMINS)  *  60  +  RLSECS)  / 
1  RHOSEC 

IF  (ILONS  .GT.  0)  3SLAL2  =  TWOPI  -  RSLAL2 

WRITE  (6,201) 

C 

C  READ  THE  RECORDS  POSITIONS  AND  OBSERVED  RATES,  THEN 

C  COMPUTE  THE  TOTAL  DISTANCE. 

1  14 


c 

DO  550  I-1,IRBC,1 

HEAD  (5,  104)  ILAT,ININ  ,RSEC,  ILON  ,  ILHIN  ,  RLSEC  , 

1  RATE1, RATE2 

104  F0RMAT(1X,I3,I3,F6.3,  15,13 ,F6. 3 , 3X,  2F1 0. 2) 

C 

C  CONVERT  LAT  AND  LONS  TO  RADIANS. 

C 

RPOST  =  ((IABS(ILAT)  *  60  ♦  IMIN)  *  60  ♦  RSEC)  / 

1  RHOSEC 

IF  (I LAT  .LT.  0)  RPOST  ■  -RPOST 

RPOSTL  =  ( (IABS  (ILON)  *  60  ♦  ILHIN)  *  60  ♦  RLSEC)  / 

1  RHOSEC 

IP  (ILON  .GT.  0)  RPOSTL  =  TWOPI  -  RPOSTL 
C 

C  COMPOTE  DISTANCES  AND  AZIMUTHS  PROM  THE  OBSERVED  POINT 

C 

C  MASTER 

C 

CALL  INVER1 (A ,F , RPOST , RPOSTL, RM AST ,R MASTL, FAZM , BAZM , 

1  DISTM) 

ODISTM  *  (RN  *  DISTM)  /  ONCOV 
CALL  SECF AC  (ODISTM, SFM) 

C 

C  SLAVE1 

C 

CALL  INVER1 (A , F , RPOST , RPO STL, RS LA V 1 , RSLAL1 , F AZ1 , BAZ 1 , 

1  DIST1) 

UDIST1  =  (RN  *  DIST1)  /  ONCOV 
CALL  SECF  AC  (ODI  ST  1 ,  S  P  1) 

C 

C  SLAVE2 

C 

CALL  INVER1 (A, F , RPOST , RPO STL,RS LAV2, RSLAL2, F AZ2 , 3 AZ2, 
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1  DIST2) 

0DIST2  =  (BN  *  DIST2)  /  ONCOV 
CALL  SECFAC  (UDIST2,SF  2) 

C 

C  DETERMINE  WHICH  LAT  AND  LONG  THE  ASP  CORHECTOB 
C  IS  ASSIGNED  TO. 

C 

CALL  ASSIGN  (RPOST ,  RPO  STL,  JN  1,  JN2) 

C 

C  COMPOTE  THE  RATES  AND  COMPARE  TO  THE  OBSERVED  RATES 
C 

TDM1  =  ODIST1  -  ODIST M  +  SF 1  -  SFM  ♦  DELAY 1 
DIFP1  =  TDM1  -  RATE1 
C 

TDM2  =  ODIST2  -  ODIST M  +  SP2  -  SFE  +  DELAY 2 
DIFF2  =  TDM 2  -  RATE2 
C 

C  SOM  THE  DIPPERENCES  TO  THE  MATRIX  AND  COONT  THE  NOMBER 

C  OF  ASF  CORRECTORS  FOR  EACH  BLOCK  TO  LATTER  DETERMINE 

C  THE  MEAN. 

C 

ASFCR1  (JN1,  JN2)  *  ASF  CR 1  ( JN  1 ,  JN  2)  ♦  DIFF1 
ASFCR2  ( JN1 ,  JN2)  *  ASFCR2  { JN  1 ,  JN2)  ♦  DIFF2 
INOI  (JN1,JN2)  *  INOI  ( JN1,  JN2)  ♦  1 
550  IN02  (JN1,  JN2)  =  IN02  ( JN1,  JN2)  ♦  1 
C 

C  DETERMINE  THE  MEAN  ASF  CORRECTOR  FOR  EACH  LAT  AND  LONG 

C 

DO  650  1=1,16,1 
WRITE  (6,201) 

DO  600  J* 1,26,1 

IF  (INOI  (I, J)  .EQ.  0)  GO  TO  580 

ASFC R 1  (I,  J)  =  ASFCR1  ( I,  J)  /  INOI  (I,  J) 

580  IF  (IN02  (I ,  J)  .EQ.  0)  GOTO  590 
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ASFCR2(I,J)  »  ASFCR2(I,J)  /  IN02(I,J) 

590  IF  (ASFCR1  (I,  J)  .SQ.  0.0)  ASFCR1(I,J)  »  9.99 

600  IF  (ASFCR2  (I,  J)  .EQ.  0.0)  ASFCR2(I,J)  *  9.99 

650  CONTINUE 
C 

C  WHITE  THE  CORRECTORS  IN  MATRIX  FORMAT 

C 

C 

IF  (I D1  .EQ.  1)  WRITE  (6, 1000) 

IF  (ID1  .EQ.  2)  WRITE  ( 6 , 100  1) 

IF  (ID1  .EQ.  3)  WRITE(6,  1002) 

WRITE  (6, 201) 

C 

WRITE  (6,  *00)  (  (ASFCR1  (I,J)  ,1-1,16)  ,J*1,26) 

DO  7  00  L-1,5,  1 
WRITE  (6,201) 

201  FORMAT  (1H  ) 

700  CONTIHOE 
C 

IF  (I D2  .EQ.  1)  WRITE  (  6,  1000) 

IF  (ID2  .EQ.  2)  WRITE  ( 6 , 100  1 ) 

IF  (ID2  .EQ.  3)  WRITE  (  6,  100 2) 

WRITE  (6,201) 

WRITE  (6,300)  (  (A  SFCR2  ( I,  J)  ,1  =  1 , 1 6)  ,  J  =  1 , 26) 

300  FORMAT  ( 16F6 .2) 

1000  FORMAT  ( '  TABLE  FOR  9940-W  ') 

1001  FORM  AT  ( '  TABLE  FOR  99  40-X  ') 

1002  FORMAT  ( '  TABLE  FOR  9940-T  *) 

STOP 

END 

SUBROUTINE  SECF  AC(  UTD  1ST,  S F) 

C - 

C 


117 


4 


C  THIS  ROUTINE  HILL  COMPOTE  THE  SEA  SECONDARY  PACTOR 

C  OTDIST  =  TOTAL  DISTANCE 

C  SF  =  SECONDARY  PACTOR 

C 

IMPLICIT  REAL*8  (A-H,  O-Z) 

C 

C 

C  COEFFICIENTS 

C 

AO  =  129.04323 

A1  =  -0.40758 
A2  *  0.00064576813 

C 

BO  =  2.741282 

B1  =  -0.011402 
B2  =  .000327748  15 
C 

IF  (UTDIST  .GT.  537)  GO  TO  10 

SF  =  (  BO  /  OTDIST)  ♦  B1  ♦  (  B2  *  OTDIST) 

GO  TO  20 

10  SF  =  (  AO  /  OTDIST)  ♦  A1  +  (  A2  *  OTDIST) 

C 

20  CONTINUE 

RETURN 
END 

SUBROUTINE  ASSI GN ( RL A T, RLON, I, J) 

C - 

c 

C  SOBROOTINE  ASSIGN  SELECTS  THE  COLOMN  AND  ROW  FOR 

C  THE  LATITUDE  AND  LONGITUDE  OF  THE  RECORD. 

C  RLAT  =  POSITION  LATITUDE  IN  SECONDS 

C  RLON  =  POSITION  LONGITUDE  IN  SECONDS 

C  I  *  COLOMN 
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IMPLICIT  REAL*8  (A-H,0-Z) 

DATA  RHOS EC/ PI/2.0  62  6  4  8  06  2  4  7 IDO 5/ 3. 1 415926535898D0/ 

TWOPI  *  PI  *  2.0 

CONVER  RLAT  A  HD  RLON  TO  SECONDS 

HLAT  =  RLAT  *  RHOS  EC 
HLON  »  TWOPI  -  RLON 
HLON  =  HLON  *  RHOS  EC 

STARTING  LAT  AND  LONG  FOR  SEARCH 

LAT  =  37/05/00.0  LONG  a  122/05/00.0 

SLAT  =  133200.0 
SLON  =  439500.0 

LATITUDE 

J  =  0 

DO  10  IC  =  1,26,1 
SLAT  =  SLAT  -  6  0.0 
RULAT  =  SLAT  ♦  30.0 

RLLAT  =  SLAT  -  30.0 

IF  (HLAT  .LT.  RULAT  .AND.  HLAT  .  GE.  RLLAT)  GO  TO  15 
J  *  IC 

LONGITUDE 

1*0 

DO  20  IC  =  1,16,1 
SLON  *  SLON  -  6  0.0 
RLLON  *  SLON  ♦  30.0 


SUBROUTINE  INVER1 (A,  ? INV, GLAT1 , GLON 1 ,GLAT2 , GLON2 ,FAZ 
1  B  AZ  ,  S) 


See  Appendix  B  for  subroutine  INVER  1 


DATA  SET  EXAMPLE  -  PROGRAM  ASFSEL 


TABLE  FOR  9940-Y 

50* 


1220  /  00* 


[i 


-0.41 

-0.53  -0.44 

45*  -0.56  -0.51  -0.46  -0.40 

-0.39  -0.54  -0.39  -0.55  -0.45 


0.63  -0.57  -0.49 
-0.54  -0.50 


1210  /  50* 


TABLE  FOR  9940-Y 
121®  /  55 « 


-0.45  -0.47  -0.41  -0.37  -0.60 


-0.51  -0.50 


-0.43  -0.44  -0.47  -0.44 


-0.38  -0.49  -0.50  -0.46  -0.52  -0.46 
-0.54  -0.51  -0.49 
-0.54  -0.57  -0.49 
-0.55  -0.55  -0.50 
-0.67  -0.61  -0.50 
-0.64  -0.65 
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